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4.  Introduction 

At  the  present  time  the  best  way  to  fight  breast  cancer  is  through  early  detection.  The 
size  of  malignancies  in  the  early  stage  is  very  small,  and  the  absorption  of  x  rays  by  the 
malignancies  and  the  breast  tissue  is  similar.  Thus  high  contrast  and  good  resolution  are 
especially  important  in  mammography.  In  current  digital  mammography  systems,  contrast  loss 
occurs  from  scattered  x  rays  that  strike  the  detector  from  all  directions,  and  resolution  is  limited 
by  the  pixel  size  of  the  digital  detector  for  computed  radiography  phosphor  plates  (CR  plates), 
or  the  conversion  efficiency  and  total  number  of  pixels  for  visible  light  CCD  systems. 

Capillary  x-ray  optics  use  total  external  reflection  to  guide  x  rays  through  tiny  hollow 
glass  tubes,  similar  to  the  way  conventional  fiber  optics  transmit  light.  A  capillary  lens  is  a 
bundle  of  hollow  glass  tubes  with  channel  sizes  as  small  as  few  microns.  The  critical  angle  for 
total  external  reflection  is  as  small  as  a  few  mrad  (1.5mrad  for  20  keV  photons).1,2  A  capillary 
lens  can  provide  almost  total  scatter  rejection.  A  well-designed  lens  can  also  provide 
magnification  without  image  degradation  from  a  finite  source  spot  size.  Thus  the  effective 
resolution  can  be  improved  while  a  CR  plate  is  used.  Capillary  optics  can  also  demagnify  the 
image  to  match  a  direct  area  detector,  such  as  a  CZT  detector3,4,  which  maybe  available  in  the 

near  future.  Unlike  the  fiber-optic-CCD  combination,  this 
system  does  not  involve  the  use  of  a  phosphor  to 
convert  the  x  rays  to  visual  light.  The  solid  state  detector 
can  provide  almost  100%  quantum  efficiency. 

Capillary  optical  mammography  has  been 
cooperatively  investigated  by  Center  for  X-ray  Optics  at 
University  at  Albany  and  the  Radiology  Dept,  of 
University  of  Wisconsin  for  many  years.  Previous  work 
has  demonstrated  the  bright  future  of  capillary  optics  in 
digital  mammography  with  both  contrast  and  resolution 
improvement.  Among  the  problems  found  in  the 
previous  work  are  the  small  size  of  the  lens  and  non- 
Figure  1  a)  Image  hluring  in  uniform  transmission.  Artifacts  caused  by  non-uniformity 

normal  airgap  magnification,  b)  0f  transmission  were  observed  in  the  image.  Eliminating 

Capillary  lens  between  the  patient  artifacts  by  image  processing  may  be  a  better  and  more 

and  detector  eliminates  the  focus  economic  way  than  fabricating  an  ideal  lens.  Low 

spot  hluring.  reproducibility  can  be  a  problem  in  lens  manufacture.  It 

is  extremely  important  to  design  a  large  lens  by 
combining  several  small  pieces.  So  a  better  understanding  of  lens  behavior  is  necessary. 
Theoretical  simulation  has  been  proved  to  be  very  helpful  in  understanding  the  capillary  and 
lens  behavior,  and  provide  feedback  to  the  manufacturing  process.  It  would  also  be  easier  to 
manufacture  large  lenses  if  they  were  shorter;  for  that  reason,  lead  glass  optics,  which  would 
stop  scattered  radiation  with  a  shorter  length,  were  investigated. 

The  proposed  work  is  focused  on  lens  measurements,  computer  simulations  and  image 
processing.  During  the  first  year  of  the  work,  a  lens  simulation  has  been  developed  and  applied 
to  the  experimental  data.  The  simulation  has  proved  to  be  helpful  to  understand  the  lens 
behavior.  MTF  is  usually  measured  by  taking  an  image  of  a  slit.  It  is  difficult  to  measure  a 
stationary  lens  with  this  method  when  the  transmission  of  the  lens  is  not  uniform.  A  method  to 
measure  MTF  by  taking  an  edge  image  was  developed  with  background  deduction  technique. 
The  lens  simulation  was  expanded  in  the  second  year,  and  lead  glass  optics  were  investigated. 
During  the  second  year,  the  principal  investigator  was  changed  from  Dr.  Lei  Wang  to  Dr.  Hui 
Wang,  and  the  third  year  the  principal  investigator  was  changed  from  Dr.  Hui  Huang  to  Dr.  Suparmi. 
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Figure  2.  Experimental  setup. 


Figure  3  Electronic  system. 


5.  Body 

5.1.  DEFECT  SIMULATION 

An  important  step  in  optics  development  was  the  breakthrough  in  modeling 
capacity,  which  has  led  to  an  unprecedented  level  of  theoretical  understanding  of  the 
basic  properties  of  polycapillary  fibers,  especially  fiber  defects.  To  evaluate  the 
experimental  performance  of  polycapillary  fibers,  and  design  capillary  optics,  it  is 
necessary  to  be  able  to  predict  theoretical  behavior  for  complex  geometries.  The 
modeling  program  for  single  fibers  is  based  on  a  Monte  Carlo  simulation  of  simple 
geometrical  optics.  The  computational  speed  is  greatly  enhanced  by  a  reduction  to  two 
dimensions  by  projecting  the  trajectory  onto  the  local  fiber  cross-section.7  Reflectivities 
are  computed  from  standard  tables.8  Significant  recent  progress  has  been  made  in 
understanding  the  effect  of  capillary  profile  error,  waviness,  and  roughness  on  the 
transmission  spectra.9  This  is  extremely  important  in  providing  feedback  to  the 
manufacturing  process.  It  has  been  found  that  extremely  good  fits  can  be  produced  with 
only  two  fitting  parameters  if  a  more  physical  model  of  waviness  is  employed.10 

The  experimental  arrangement  for  single  capillary  measurements  is  shown  in 
Figure  2  and  Figure  3.  An  optical  rail  affixed  to  an  optical  table  carries  an  x-ray  source, 
fiber  platform,  and  x-ray  detector.  Each  can  be  positioned  independently  in  three 
dimensions.  A  collimator  is  placed  before  the  fiber  and  any  remaining  x-ray  leakage 


Type 

Description 

Outer  Diameter, 
mm 

Channel  Size, 
pm 

Open  Area 

Length, 

mm 

■ 

Borosilicate 

0.5 

12 

65% 

105 

B 

Lead  glass 

0.5 

12 

52% 

95 

C 

Borosilicate 

0.75 

22 

50% 

136 

D 

Borosilicate 

4 

12 

55% 

130 

E 

Borosilicate 

0.3 

4-5 

55% 

105 

Table  1.  Description  of  polycapillary  fibers. 
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■  exp,  Fiber  a  . sim,  Fiber  a  around  the  fiber  is  eliminated  with 

metal  powder  or  filings.  The 
measured  fibers  are  described  in 
Table  1. 

The  results  of  transmission  studies 
as  a  function  of  photon  energy  are 
shown  in  Figure  4. 11,12  All  of  the 
fibers  except  the  lead  glass  have 
transmissions  at  20  keV  nearly  equal 
to  their  fractional  open  area  (the 
fraction  of  the  cross  section  of  the 
capillary  which  is  open  space,  the  rest 
being  glass  walls).  This  transmission 
corresponds  to  the  primary 
Photon  energy  (KeV)  transmission  expected  for  a  linear 

Figure  4.  Measured  transmission  versus  energy  for  Capillary  Optic  employed  as  an 

polycapillary  fibers  listed  in  Table  1.  antiscatter  grid. 

The  transmission  falls  off  at  higher 
energies.  Modeling  this  transmission  spectra  has  lead  to  a  greater  understanding  of 
defects  in  polycapillary  optics. 

In  Figure  5,  simulations  with  or  without  roughness  corrections  are  compared  with 
the  experimental  data.  In  these  measurements  the  source  is  scanned  transverse  to  the 
fiber  axis.  The  simulation  with  a  roughness  height  of  0.5  nm  fits  the  experimental  data 
quite  well.  It  is  definitely  over-corrected  when  the  roughness  is  1 .0  nm.  The  same 
simulations,  shown  in  Figure  7,  are  also  carried  out  at  68  keV,  where  the  width  of  the 
curve  is  narrower  than  that  at  10  keV  because  of  the  smaller  critical  angle.  As  a  result, 
photons  also  experience  fewer  bounces  on  average.  Simulations  with  roughnesses  as 


•  data 


Source  position(mm) 

Figure  5.  Source  scans  at  10  keV. 


Figure  6.  Transmission  spectra  of  Fiber  3  simulated 
with  different  bending  curvature  alone  and  compared 
with  experimental  data 
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■  data 

_ z=0 

.  z=1  nm 


Source  position(mm) 


Figure  7.  Source  scans 
at  68  ke  V. 


large  as  1.0  and 
2.0  nm  still  could 
not  fit  the  data. 

Knowing  that  1.0 
or  2.0  nm 
roughness  is 
definitely  too 
large  at  10  keV, 
we  can 

determine  that 
the  roughness 
correction  by 
itself  is  not 
sufficient  at  high 
energy  to 
reproduce 
source  scan 
measurements. 

Other  effects 
need  to  be 

considered.  These  are  bending  and  waviness. 


Figure  8.  Simulations  of  transmission  spectra  with 
waviness  only  compared  with  the  experimental  data. 


Figure  9.  Simulated  source  scan  curves  compared  with  experimental  data  at  four  different  photon  energies. 
Parameters  are:  R  =  125  m,  A9max=  0.35  mrad,  roughness  height=0.5  nm. 
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A  slight  bending  of  the 

capillary  can  dramatically  reduce  50 .  .  .  m  .  . 

the  transmission  of  high  energy  *  "  •  • . 

photons  because  of  the  small  ^  . 

critical  angle.  A  comparison  ’  v  • 

between  experimental  data  and  ^  N  • 

simulations  with  different  ^  ^0  -  s"V . 

bending  is  shown  in  Figure  6.  §  '"V 

The  figure  shows  that  the  |  35  -  ■ 

simulations  with  bending  alone  £  .  ■  exp  data,  Fiber  C  \ 

do  not  fit  the  experimental  data  §  30 - ami,  R=125rr\  M  rrax=0.35macl  \ 

well  which  indicates  that  £  . sintf,  R=110m,  Ae  rax  =0.25nrad 

bending  is  not  the  only  factor 

which  causes  the  high  energy  ~  ”  R^Om,  A0  max  =0.45mrad 

transmission  to  drop.  However,  ' 

from  Figure  6,  we  can  see  that  20q - 1 - 20 - * - 40 - ’ - 60 - * - j 

the  range  of  the  bending  radius  ! 

must  be  larger  than  1 00m  to  give  Phcton  energy(Ke\/) 

the  observed  transmission  at  the 

highest  energy(80KeV)  Figure  10.  Simulated  transmission  spectra  with  different 

Capillary  surface  bending  and  waviness  compared  with  the  experimental  data  in 

oscillations  with  wavelengths  search  for, he  bes,  fining  of  Fiber  C. 

shorter  than  the  capillary  length 

and  longer  than  the  wavelength  of  the  roughness  are  called  waviness.  The  detailed 
shape  of  waviness  is  unknown.  Its  average  effect  can  be  considered  as  a  random  tilt  of 
the  glass  wall,  so  that  the  grazing  angle  of  the  photon  is  changed  by  a  random  amount 
86,  after  every  bounce.  $9  is  a  random  number  between  -A9max  and  A0max  if  0>zl<9max. 
The  maximum  random  tilt  angle  A0max  is  an  adjustable  parameter  which  depends  on 
the  waviness  of  the  polycapillary  fiber  To  keep  ^'positive,  56  is  taken  to  be  a  random 
number  between  -6  and  A0max  when  9  <A6max .  Since  a  photon  with  a  incident  angle 
smaller  than  zl<9max  has  a  larger  chance  to  experience  an  angle  increase  than  an  angle 
decrease,  this  is  physically  reasonable.  In  Figure  8,  simulations  with  waviness 


■  exp  data,  Fiber  C 

—  siml,  R=125m,  ag  max=0.35mrad 
•  •  •  si  m2,  R=110m,  ABmax  =0.25mrad 

-  sim3,  R=170m,  ag  max  =0.45mrad 


Photon  energy(Ke\/) 

Figure  10.  Simulated  transmission  spectra  with  different 
bending  and  waviness  compared  with  the  experimental  data  in 
search  for  the  best  fitting  of  Fiber  C. 


Figure  11.  Scheme  of  an  x-ray  beam  interacts  with  a 
random  tilt  surface. 


corrections  with  A0max  set  at  1 
mrad  and  2  mrad,  which  is 
comparable  to  the  critical 
angle,  are  compared  with  the 
experimental  data.  This  figure 
shows  that  simulations  with 
waviness  alone  do  not  fit  the 
experimental  data.  This  is 
because  the  waviness 
correction  changes  the 
reflected  angle,  not  the  profile. 
In  fact  the  capillary  is  still 
considered  to  be  straight,  so 
those  photons  which  have  few 
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reflections  will  not  be  significantly 
effected  by  waviness. 

Finally  the  waviness  and 
bending  are  combined  by  increasing 
the  bending  radius  R,  roughly 
determined  in  Figure  6,  and  adding  a 
waviness  parameter,  ABmax.  Several 
trials  are  shown  in  Figure  10.  Sim2  has 
too  much  bending  and  not  enough 
waviness;  sim3  has  too  much  waviness 
and  not  enough  bending;  siml  is  the 
best  fit.  Roughness  is  also  included  in 
those  simulations.  The  source  scan 
simulation  with  the  three  fixed 
parameters  are  plotted  along  with  the 
experimental  data  in  Figure  9  for  four 
more  photon  energies.  They  all  fit  quite  well.  This  model,  hereafter  labeled  Ml, 
assumed  a  uniform  distribution  of  tilt  angles,  which  is  not  quite  physical.  In  more  recent 
work,  discussed  as  Model  M2,  it  is  assumed  that  these  tilt  angles,  <|>,  are  normally 
distributed  in  the  range  (-ti/2,  7r/2)  with  the  mean  value  equal  to  zero.  For  high  quality 
optics,  the  standard  deviation  of  this  normal  distribution,  a,  is  much  smaller  than  the 
critical  angle,  0C.  The  probability  distribution  of  tilt  angles,  <j>,  is 


Figure  12.  Three  surfaces,  OA1,  OA2  and  OA3,  with 
different  tilt  angles,  <j>l,  (j>2  and  <j)3,  respectively,  from  the 
nominal  surface  00 


G(</>)  = 


(1) 


In  this  work,  Model  M2,  consideration  was  taken  of  the  fact  that  the  surface  tilt  angle  will 
affect  the  probability  of  x-ray  impact  on  that  surface.  Taken  to  extremes,  a  surface 
region  perpendicular  to  the  beam  is  much  more  likely  to  intercept  the  beam  than  a 
surface  region  parallel  to  the  beam.  Figure  12  displays  three  surfaces,  OAi,  OA2  and 
OA3,  with  different  tilt  angles,  <j>i,  <j)2  and  <j>3,  respectively,  from  the  nominal  surface  OO’. 
The  projections  onto  the  nominal  surface  for  the  three  surfaces  are  equal,  OA^cos^i  = 
OA2*cos<)>2  =  OA3*cos<J)3  =  OO’.  The  probability  of  incidence  of  a  parallel  x-ray  beam  with 
incident  angle  0  (with  respect  to  the  nominal  surface  OO’)  hitting  these  tilted  surfaces  is 
given  by  their  corresponding  perpendicular  length  AjBj  (for  j  =  1,  2  or  3) 

P,  cc  A ,B  =  OA,  -  sin (0  +  ^)=  ----  •  sin(i9  +  (j), ),  ^ 

COS  <j)j 

We  call  this  the  tilt-corrected  probability  distribution.  The  complete  description  of  this 
distribution  is 


F  sin(  0  +  <f>) 
cos(  </>) 


-6  <<j>  <  — 
2 

2 


H  (<!>)  =  { 


0, 


(3) 
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where  0  is  the  incident  angle,  <j)  is  the  tilt  angle  and  F  is  a  normalization  constant. 

Combining  the  normal  distribution  G(<j>)  with  the  tilt  correction  H(<j))  gives  (j>  for  a 
certain  incident  angle  0  as  follows: 


sin(<9+(Z>) 

cos^ 

0, 


-^<f, 
f ^ 


where  0  is  the  incident  angle,  <j>  is  the  tilt  angle  and  K  is  a  normalization  constant. 
Noting  that 

sin^+^)  .  _  ,  (5] 

- —  =sln0^-cos6,•  tan^ 

cos <j) 

i 

and  that  both  the  incident  angle,  0,  and  the  tilt  angle,  (j>,  are  very  much  less  than  40 
mrad,  we  use  the  approximation 

f  z*L 

P(t)  =  \  2, 

0, 


(4) 


(6) 


where  the  normalization  constant ,  C,  is 

C  «  —7= — - - 

'JIkOct  +  2a 2 

to  speed  the  calculation.  P(<j>)  is  called  the  tilt-corrected  normal  distribution.  In  model 
M2,  c  is  the  variable  that  describes  the  amount  of  waviness.  An  example  of  J((|))  and 
P(<j>)  is  shown  in  Figure  13.  The  result  of  the  waviness  calculation  is  shown  in  Figure  14 
The  model  M2  shows  a  much  larger  drop  at  high  photon  energies  for  lower  values  of 
the  waviness  than  does  model  Ml. 
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Figure  13.  P((j))  andJ(<f>)  versus  tilt  angle  <t> 
when  incident  angle  6  =  0. 009rad, 
standard  deviation  cr=  0. 0004rad.  tj>  is  in 
rad. 


A  comparison  of  the  simulation  with  bending  and 
waviness  to  the  experimental  data  is  shown  in 
Figure  15. 

Using  a  similar  technique,  the  simulation 
was  compared  to  the  fibers  listed  in  Table  1, 
using  the  parameters  listed  in  Table  2.  The 
results  are  plotted  in  Figure  16  and  Figure  17.  In 
Figure  16,  the  transmission  for  fiber  A  shows  a 
rapid  drop  for  energies  above  30  keV.  Although 
fiber  D  has  lower  fractional  open  area  than  fiber 
A,  its  transmission  exceeds  that  of  fiber  A  at 
energies  above  30  keV.  This  is  because  fiber  A 
is  thin  (0.5  mm  in  outer  diameter)  and  flexible, 
therefore  difficult  to  keep  straight  in  the 
measurement  apparatus.  Model  M2,  which  can 
vary  only  bending  and  waviness,  requires  a  much 
sharper  bend  for  fiber  A  than  for  fiber  D.  This 
bending  is  more  significant  at  high  energies, 
where  the  critical  angles  are  smaller. 
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Figure  14.  Simulations  of 
transmission  spectra  for  fiber  C  with 
only  waviness  compared  with  the 
experimental  data.  This  figure  shows 
the  effects  of  waviness.  It  also  shows 
simulations  using  waviness  alone  do 
not  fit  the  data.  The  simulations  do 
not  include  roughness  or  bending. 


Figure  15.  Simulated  transmission  spectra  using 
model  M2  with  fixed  waviness  (cr=0.2  mrad)  and 
different  bending  (bisection  method),  compared 
with  the  experimental  data.  The  simulation  curves 
shown  are  just  a  few  representative  cases  in  the 
bisection  process.  For  M2,  the  best-fit  parameters 
are  <j=  0.2  mrad  and  R  =  225  m;  for  Ml,  they  are  z 
=  0.5  nm,  s  =  6  pm,  A0max  =  0.35  mrad  and  R  =  125 
m. 
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Transmission  curves  for  fiber  C  and  D 
are  similar  in  shape  and  are  nearly  flat 
up  to  60  keV,  as  shown  in  Figure  17. 
However,  the  smaller  channel  size  for 
fiber  D  results  in  the  simulation  being 
less  sensitive  to  bending.  Therefore  a 
larger  bending  curvature  (smaller 
radius)  is  required  for  fiber  D.  The 
waviness  correction  for  fiber  D  must 
therefore  be  smaller  to  keep  the 
transmission  about  the  same  as  for 
fiber  C. 


Photon  energy  (KeV) 

Figure  16.  Simulations  of  transmission  spectra  of fiber 
A  and  fiber  D  with  their  best-fit  parameters  compared 
with  the  experimental  data.  All  parameters  related  to 
this  figure  are  listed  in  Table  1. 


Fiber  E  is  the  thinnest  fiber. 
Table  2  shows  that  it  has  the  largest 
bending  correction  because  of  its 
flexibility.  Its  transmission  curve  is  flat 
up  to  40  keV.  This  is  due  to  its  very 
small  channel  size,  as  discussed  in  the 

comparison  of  fiber  C  and  fiber  D.  However,  if  the  channel  size  is  too  small,  as  in  fiber 
E,  it  also  results  in  more  reflections  being  needed  for  a  photon  to  traverse  the  fiber  and 
may  have  introduced  other  defects  such  as  blocked  channels.  This  is  why  the 
transmission  is  only  40%  for  the  energies  below  40  keV  although  the  open  area  is 
around  55%. 


The  new  model  yields  as  high  a  quality  fit  with  two  fewer  fitting  parameters  than  used  in 

model  Ml  in  Figure  4.  The  high 
transmission  and  the  simulation 
results  show  that  the  quality  of  the 
capillary  fibers  is  quite  good.  The 
bending  radius  is  above  130 
meters  for  type  3  and  type  4 
capillaries.  It  is  hopeful  that  we 
can  further  improve  the  high- 
energy  transmission  performance 
of  polycapillary  fibers  by 
decreasing  to  an  optimum 
channel  size,  and  making  them 
more  rigid. 


Photon  energy  (KeV) 

Figure  1 7.  Simulations  of  transmission  spectra  of  fiber  C, 
fiber  D  and  fiber  E  with  their  best-fit  parameters  compared 
with  the  experimental  data.  All  parameters  related  to  this 
figure  are  listed  in  Table  2. 


tranmission10la#1 


Simulations,  Measurements  and  Image  Processing  for  Capillary  Optical  Digital  Mammography 


13 


Fiber  Description 

Model  Ml 

Model  M2 

Fiber 

# 

Type 

Outer 

Diameter 

(mm) 

■m 

Area 

Length 

(mm) 

z 

nm 

s 

fim 

R 

m 

A0max 

mrad 

R 

m 

c 

Mrad 

A 

1 

0.5 

12 

65% 

105 

0.7 

6 

105 

0.4 

60 

0.225 

C 

3 

0.75 

22 

50% 

136 

0.5 

6 

125 

0.35 

225 

0.2 

D 

4 

4 

12 

55% 

130 

0.8 

6 

110 

0.285 

139 

0.125 

E 

5 

0.3 

MEM 

55% 

105 

0.7 

6 

28 

0.2 

31 

0.09 

F 

4 

4 

12 

55% 

130 

0.8 

6 

90 

0.45 

90 

0.18 

Table  2  Parameters  for  best-fit  simulations.  R  is  the  bending  radius.  For  Ml,  AOmaxis  the 
amount  of  waviness,  z  is  the  roughness  height  and  s  is  the  roughness  correlation  length. 
For  M2,  <7  is  the  standard  deviation  of  the  waviness. 


Fiber  B  in  Figure  4  has  poor  transmission;  it  is  a  lead  glass  fiber,  with  high 
surface  roughness  and  waviness.  More  recently,  higher  quality  lead  glass  fibers  have 
been  produced.  The  transmission  as  a  function  of  photon  energy  for  these  fibers  is 
shown  in  Figure  18.  These  leaded  fibers  have  higher  absorption  of  high  angle  photons 
than  the  lower  density  borosilicate  glass  fibers,  and  therefore  would  provide  even  better 
scatter  rejection,  as  is  described  in  section  0. 


Transmission  was  also  measured  of  deliberately  curved  fibers,  which  can  be 
used  to  deflect  x-ray  beams  over  larger 
angles.  The  results  are  shown  in 

Figure  19.  x  _ 
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Figure  19.  Transmission  efficiency  of 350  /mu  OD 
polycapillary  fiber  with  6  pm  channel  diameter.  The 
fiber  length  is  25  cm.  Simulated  values  are  shown  for 
20,  30,  40,  50  keVxrays.  The  parameters  are 
curvature,  R;  waviness,  w;  and  roughness,  z. 


Figure  18.  Recently  manufactured  lead  glass 
fibers. 
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Figure  20.  Measured  local  divergence  of  the  output  of  multi-fiber  I  at  8 
keV. 


5.2.  SIMULATION  OF 
MULTIFIBER  OPTICS 

Measurements  of 
the  exit  divergence  of  the 
collimating  optic  have 
been  performed  at  8  keV 
and  at  20  keV.  Exit 
divergence  is  important  if 
a  capillary  optic  is  to  be 
used  as  a  collector  fore 
slit  before  the  patient.  In 
this  case,  the  spatial 
resolution  will  depend  on 
the  angular  spread  of 
x-rays  at  the  exit  end  of 
the  polycapillary  fibers. 

Large  angular  divergence  is  equivalent  to  a  large  focal  spot  size  in  a  conventional 
system. 

The  divergence  can  also  be  used  to  assess  fiber  alignment  in  the  lens.  The 
beam  exiting  the  optic  is  characterized  by  a  divergence  that  arises  from  two  effects.  The 
first  is  the  local  divergence  of  the  x-rays  emerging  from  each  channel,  usually  between 
0C,  the  critical  angle  for  total  reflection  and  20c.  There  is  also  fiber  misalignment,  the 
deviation  of  the  individual  channel  axis  direction  from  the  optic  axis  direction.  Local 
divergence  depends  upon  the  surface  roughness,13  waviness,  and  bending  of  fibers 
through  the  optic.  Figure  20  shows  the  measured  local  divergence  of  the  output  of  the 
collimating  lens  at  8  keV  with  the  aperture  at -15  cm,  -10  cm,  0  cm,  10  cm  and  15  cm 
off  the  axis  of  the  lens.  The  FWHM  of  each  divergence  curve  is  in  the  range  of  3.8  mrad 
to  3.9  mrad,  which  is  very  close  to  the  critical  angle  of  about  4  mrad  at  8  kev.  The  peak 
centers  of  these  five  local  divergence  curves  are  systematically  shifted.  This  shift  could 
be  caused  by  the  output  ends  of  the  fibers  inside  the  lens  being  slightly  convergent 
rather  than  parallel.  The  widths  of  the  measured  local  divergence  curves  shown  in 

Figure  20  are  larger  than  that 
predicted  by  an  ideal  lens 
divergence  simulation,  shown  in 
Figure  21. 

The  divergence  of  the  modeled 
ideal  lens  is  low  because  the 
nearly  straight  central  fibers,  if 
ideal,  would  not  increase  the 
divergence  above  the  entrance 
divergence  due  to  the  source 
spot  size.  Figure  23  shows  the 
simulated  divergence  profile  of  x- 
Divergence  angle  (mrad)  rayS  exjtjng  from  an  ideal  straight 

Figure  21.  An  ideal  whole  lens  simulation  at  8  keV gives  a  fiber,  which  has  a  FWHM  of  2.36 

divergence  FWHM  of  2. 5  mrad.  m  rad . 
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Angle  (mrad)  Angle  (mrad) 


Figure  23.  For  an  ideal  fiber  (with  no  surface  Figure  22.  A  simulation  with  a  waviness  of  0.15 

defects),  the  divergence  simulation  yields  a  mrad  increases  the  divergence  value  from  2.4 

FWHM  of  2.36  mrad  at  8  ke  V.  mrad  for  an  ideal  straight  fiber  to  4. 0  mrad. 

The  simulations  in  Figure  21  and  Figure  23  did  not  include  the  effects  of  profile 
defects  and  waviness.  Waviness  will  increase  the  angle  of  reflection  for  x-ray  photons 
for  most  bounces  inside  the  channel.  Consequently  the  divergence  from  the  lens 
increases.  Figure  22  shows  a  simulated  transmission  of  a  straight  fiber  as  a  function  of 
exiting  angle,  using  a  waviness  of  co  =  0.15  mrad.  The  width  of  the  simulated  divergence 
curve  for  this  straight  fiber  is  about  4.0  mrad,  which  is  close  to  the  local  divergence  at 
the  center  of  the  lens  shown  in  figure  16.  Waviness  of  the  channels  changes  the 
divergence  of  the  x  rays  exiting  from  the  center  of  the  lens  from  about  2.4  mrad  to  4.0 
mrad. 


5.3.  MAGNIFYING  LINEAR  MONOLITHIC 
ANTISCATTER  OPTICS 
5.3.1.  Transmission 

Linear  monolithic  optics  taper  as  a  single  glass 
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Figure  25.  Transmission  of  tapered  lens. 
Solid  line  is  simulation. 


Figure  24.  2-D  source  scan  of  early  prototype 
optic  (input  diameter  =  2  mm,  output  diameter 
=  3  mm,  transmission  =  5  %  at  20  keV,  focal 
distance  22  cm)  This  lens  transmits  poorly 
because  only  the  central  capillary  channels 
transmit  while  the  outer  capillaries  do  not 
transmit  at  all.  Only  about  33  %  of  the  optic 
cross  section  transmits  at  5  %  or  greater. 
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unit  from  input  to  output  as  shown  in 
Figure  28.  These  optics  are  relatively 
difficult  to  manufacture.  Many  early 
prototypes  have  poor  transmission,  as 
shown  in  Figure  24.  To  help  understand 
the  defect  problem  an  extensive 
simulation  and  measurement  study  was 
performed  on  several  of  the  prototypes. 

Transmission  of  a  monolithic 
taper  is  plotted  versus  source  lens 
distance  in  Figure  25  along  with  the 
simulation  values.  The  maximum  point 
for  transmission  is  around  450  mm. 

The  second  method  can  be  affected  by 
the  source  instability. 

Transmission  was  measured  in 
a  wide  energy  range  for  the  whole  lens 
as  well  as  the  central  part  of  the  lens. 
The  transmission  spectrum  for  the 
whole  lens  is  plotted  in  Figure  26.  The 
figure  shows  that  the  transmission 
drops  quickly  for  photons  with  energy 
higher  than  20  keV,  but  the 
transmission  under  20  keV  is  almost 
flat.  The  low  transmission  in  the  high 


Photon  energy  (keV) 

Figure  26.  Transmission  of  whole  tapered  lens. 


•  Pinhole  size:  1 
mm. 

°  Pinhole  size:  0.5 
mm. 

—Simulated 
transmission, 

1  mm,  no  waviness 

0  20  40  60  80  100 
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Figure  27.  Transmission  of  the  central  part  of  the 
tapered  optic. 


energy  range  can  be  explained  by  the  non-linearity  of  the  lens.  The  outermost  channels 


are  more  bent  than  the  central  channels,  so  they  have  low  or  zero  transmission  at  the 


higher  photon 
energies.  This  is 
confirmed  by  the 
higher  transmission 
measurement  at  the 
central  part  of  the 
lens,  shown  in 
Figure  27  and  the 
transmission 
uniformity 


Figure  28.  Small  optic  scanning  geometry  for  contrast  and  imaging 
measurements. 


measurement  at 
different  energies, 
shown  in  Figure  29. 


Transmission 


uniformity  of  the  lens  was  measured  by  putting  a  200  pm  pinhole  in  front  of  the  lens  and 
performing  a  two  dimensional  scan.  At  8  keV,  the  whole  lens  transmits,  although  we 
can  see  the  transmission  is  not  completely  flat  due  to  the  non-linearity  and  the  defects 
of  the  lens,  such  as  blocked  channels.  At  25  keV,  transmission  of  the  whole  lens  drops 
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Figure  29.  Two  dimensional  scans  of  the  output  of  the  tapered  optic. 

to  25%.  Most  of  the  lens  is  still  transmitting  at  25  keV,  but  the  outer  channels  have 
lower  transmission  than  the  central  part.  At  35  keV,  the  lens  has  a  narrower 
transmitting  region.  The  measurements  show  that  the  lens  has  undesirable 
nonlinearity. 


5.3.2.  Scatter  Fraction  and  Contrast  Enhancement  Measurements 


Figure  30.  Experimental  setup  for  scatter  rejection  and  contrast 
measurements 


The  experimental  setup  for  scatter  rejection  and  contrast  measurements  is  shown  in 
figure  30.  The  x-ray  generator  used  in  the  experiment  was  a  Mo  source  with  a  300  pm 
source  spot  size  and  maximum  operating  voltage  of  35  kVp.  The  source  head  was 
mounted  on  two  stages  which  could  be  translated  in  the  two  directions  transverse  to  the 
x-ray  beam.  Thick  lead  shielding  was  used  to  prevent  background  x  rays  from  reaching 
the  imaging  plate.  The  borosilicate  monolithic  optic  was  put  on  aluminum  plate  and 
covered  by  powder  filling  to  prevent  x-ray  leakage  around  the  optic.  A  Fuji  imaging  plate 
with  50  pm  resolution  was  used  to  record  the  image.  The  Lucite  phantom  used  in  the 
contrast  measurement  was  45  mm  thick.  Several  lead  strips  with  different  width  were 
attached  to  the  phantom  to  measure  scatter  transmission. 

A  prototype  borosilicate  monolithic  optic,  166  mm  in  length,  with  input  diameter  of 
4.5  mm  and  7.6  mm  output,  was  used  to  measure  scatter  fraction  and  contrast 
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enhancement  of  a  Lucite  phantom  with  a  thickness  of  45  mm. 
whole  optic  was  35%  at  20  keV 
central  part  of  the  optic,  4  mm  in 


Figure  31.  Output  uniformity  of  a 
long  monolithic  optic,  L=166 
mm,  Dmput  =  4.5  mm,  Doulpul  = 7.6 


Figure  32.  Sketch 
paths  for  scatter 
measurement. 


mm  optic  provides  magnification, 

resolution,  as  well  as  reducing  scatter. 


The  transmission  of  the 
and  7%  at  60  keV.  A 
diameter  at  the 
output,  has  good 
transmission  up  to 
energy  40  keV. 

The  output 
uniformity  of  the  optic 
is  shown  in  figure  31. 
The  channels  of  the 
outer  part  of  the  optic, 
which  does  not 
transmit  well,  may  be 
blocked  by  glass 
inclusions.  Such  an 


and  thus  enhanced 


5.3. 2.1  Scatter  Fraction 

To  measure  the  scatter  transmission  of  the  optic,  several  lead  strips  with  different 
width  but  the  same  thickness  were  attached  to  the  phantom,  as  shown  in  figure  32.  The 
scatter  fraction,  SF,  was  calculated  using  formula 


SF  -  Is  /  ( Is  +  Ip  ) 


(8) 


where  Is  is  the  intensity  behind  the  lead  strip  (scatter  intensity)  and  lP  is  the  primary 

intensity,  which  is 
the  intensity 
around  the  lead 
strip  subtracted 
by  scatter 
intensity.  The 
optic  reduced  the 
scatter  fraction  by 
nearly  a  factor  of 
10,  as  shown  in 
figures  33a  and 
b. 


Figure  33a.  Scatter  fraction  with 
optic. 


Figure  33b.  Scatter  fraction  without 
optic. 


5. 3.2. 2  Contrast 
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Figure  34.  Sketch  of  Lucite  phantom 


The  Lucite  phantom  used  in  the  contrast 
measurements  was  45  mm  thick  and  had  four 
holes  with  the  same  diameter  of  1 .5  mm  and 
different  depths,  from  10  mm  to  35  mm,  as  listed 
in  Table  III.  The  diameter  of  the  holes  were 
smaller  than  the  central  area  the  optic  which 
transmits  well.  Each  image  was  taken  on  a 
single  hole.  A  sketch  of  the  Lucite  phantom  with 
contrast  was  measured  with  and  without  optic. 
The  images  with  and  without  the  optic  are 
shown  in  figures  35a  and  b.  The  contrast,  C,  is 


d  =20  mm 


d=35  mm 


C  =  In  (  T  / 1 )  (9) 


d=10  mm 


d=30  mm 


Figure  35b.  Lucite  phantom 
image  without  optic. 


where  I  is  the  intensity 
near  the  hole  and  I  is 
the  intensity  through 
the  hole.  The  contrast 
with  and  without  the 
optic,  and  also  contrast 
enhancement,  are 
listed  in  Table  3.  The 
contrast  enhancement 

is  the  ratio  of  the  contrast  with  and  without  the 
optic.  The  results  show  that  the  optic  increases  the 
contrast  by  a  factor  of  1 .7  for  the  shallowest  hole. 
The  images  in  figures  35a  and  b  are  shown  to 
scale.  A  magnification  of  1 .8  was  obtained  with  the 
optic.  This  long  monolithic  optic  with  small  area 
was  used  to  magnify  the  image,  reduce  the  scatter  transmission  and  enhance  the 

contrast. 


•  • 


Figure  35a.  Lucite  phantom  image 
with  optic. 


Hole 
depth,  d 
(mm) 

Contrast 

Ratio 

With 

optic 

No  optic 

35 

3.1  ±0.1 

2.4  ±0.03 

1.3  ±0.1 

30 

2.6  ±0.1 

1.9 ±0.04 

1.4  ±0.1 

20 

2.1  ±0.1 

1.4  ±0.1 

1.5 ±0.1 

10 

1.2  ±0.0 3 

0.7 ±0.02 

1.7 ±0.02 

Table  3.  Contrast  and  contrast 
enhancement 


5. 3. 2. 3  Disadvantages  of  Long 
Borosilicate  Optics 

Long  optics,  even  with  small  input  area, 
tend  to  have  defects  as  shown  in  figure  3.  A 
long  optic  also  leads  to  manufacturing 
complexity  and  expense.  Improvements  in 
the  technology  of  fabricating  monolithic 
optics  have  been  achieved.”4  Medical 
imaging  or  astronomy  applications, 
however,  require  large  area  optics.  For 
high-energy  applications,  borosilicate  optics 
need  to  be  quite  long.  Leaded  glass  would 
allow  the  optic  to  be  much  shorter  and  still 
good  give  scatter  rejection  because  of 
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strong  absorption  of  lead  glass.  The  shorter  optics  are  easier  to  manufacture  and  could 
be  cheaper.  In  the  following  section,  several  leaded  glass  fibers  were  measured  and 
analyzed  using  a  geometric  simulation. 

5.3.3  MTF  Calculation 


Modulation  Transfer  Function  (MTF)  is  the  most  fundamental  measurement  of  spatial 


resolution  used  in  radiology.  The  standard 

technique  is  to  image  a  slit,  determine  the  Figure  37.  Three-dimensional  edge  image  after 
line  spread  function  (LSF),  and  compute  background  subtraction. 
the  Fourier  transform.  The  MTF 

measured  for  the  mammography  system  with  capillary  optics,  compared  to  that  without  the 
capillary  optics,  will  give  the  resolution  improvement  and  also  can  be  used  to  diagnose  artifacts. 
The  structure  of  the  capillary  bundles  was  obvious.  In  this  case,  the  LSF  could  not  be 
obtained  by  imaging  a  slit  without  background  elimination.  The  background  deducted 
slit  image  is  the  slit  image  divided  by  the  image  with  capillary  only.  However,  this 
involves  image  registration.  The  information  for  background  that  can  be  used  for  image 
registration  is  too  little  in  the  slit  image.  An  alternative  method  is  calculating  the  LSF  as 
the  derivation  of  the  edge  spread  function  (ESF).  The  edge  image  left  enough 
background  for  registration  so  that  it  could  be  a  feasible  method  in  our  case. 

The  structure  of  the  capillary  bundles  was  obvious.  In  this  case,  the  LSF  could 
not  be  obtained  by  imaging  a  slit  without  background  elimination.  The  background 
deducted  slit  image  is  the  slit  image  divided  by  the  image  with  capillary  only.  However, 
this  involves  image  registration.  The  information  for  background  that  can  be  used  for 
image  registration  is  too  little  in  the  slit  image.  An  alternative  method  is  calculating  the 
LSF  as  the  derivation  of  the  edge  spread  function  (ESF).  The  edge  image  left  enough 
background  for  registration  so  that  it  could  be  a  feasible  method  in  our  case. 
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Figure  38.  Edge  spread  function  and  its  fitting  curve.  Figure  39.  Reprojection  of  a  two-dimensional 

edge  image  into  a  finely  sampled  ESF. 

A  registration  algorithm  (FMI-SPOMF)  proposed  by  Chen,  et  al.  14  was  used  in 
this  work.  This  is  a  method  to  match  a  two-dimensional  image  to  a  translated,  rotated 
and  scaled  image.  The  approach  consists  two  steps:  the  calculation  of  a  Fourier-Mellin 
Invariant  (FMI)  descriptor  for  each  image  to  be  matched,  and  then  matching  of  the  FMI 
descriptors.  FMI  descriptors  are  translation  invariant.  The  matching  of  the  FMI 
descriptors  is  to  find  out  the  rotation  and  scaling,  and  achieved  using  a  symmetric 


0  1  2  3  4  5  6  7 

cycles/mm 

Figure  40.  Calculated  MTF  compared  with  the 
MTFs  with  scan  optics,  and  without  optics. 


phase-only  matched  filtering  (SPOMF).  In 
our  case,  there  is  no  change  in  scaling,  but 
rotation  is  probably  involved.  When  the 
rotation  was  found,  image  translation  is 
found  by  SPOMF  method. 

Images  were  taken  with  edge  as 
shown  in  figure  36.  The  edge  was  made  of  a 
lead  plate.  The  size  of  each  image  is  fifty  by 
fifty  in  pixels.  They  were  taken  in  an 
experimental  digital  mammography  system 
with  a  computed  radiography  digital 
phosphor  plate.  The  digital  CR  plate  has 
limited  resolution,  5  Ip/mm,  and  is  usually 
considered  to  be  not  good  enough  for 
clinical  mammography.  But  the  effective 
resolution  of  the  system  could  be  improved 
by  the  magnification  of  the  capillary  optics. 
The  magnification  of  the  capillary  lens  used 
in  the  system  is  1 .86,  so  that  the  effective 
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resolution  should  be  improved  by  a  factor  of  1.86. 

The  registration  algorithm  was  implemented  in  IDL.  With  the  FMI-SPOMF  registration 
algorithm,  no  rotation  was  found.  The  translation  was  then  found  by  SPOMF  algorithm. 
The  difference  between  the  usual  correlation  method  and  SPOMF  method  is  that  it  only 
uses  the  phase  information.  The  phase-only  correlation  function  has  sharper  peak  than 
normal  correlation  function.  Sub-pixel  resolution  could  be  achieved  theoretically. 
However  it  was  found  that  the  resolution  is  not  good  enough  in  our  case.  This  may  be 
caused  by  the  relatively  small  image  size.  The  background-subtracted  result  was  further 
optimized  by  manual  shifts  in  sub-pixel  range.  The  result  after  manual  optimization  was 
shown  in  figure  37.  The  periodic  background  is  totally  gone.  A  cubic  spline 
interpolation  method  was  used  for  the  best  result  while  shifting  the  background  image 
relative  to  the  edge  image. 

The  presampling  MTF  includes  unsharpness  of  the  detector  and  the  sampling 
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Figure  41.  Schematic  of  the  experimental  setup  for  optic  characterization. 
aperture.  In  order  to  eliminate  the  aliasing,  a  finely  sampled  ESF  is  obtained  with  a 
slightly  angulated  edge  in  a  single  exposure  as  illustrated  in  figure  39. 15,16  The  angle  in 
our  edge  image  in  figure  37  was  calculated  to  be  around  6  °.  The  resultant  ESF  is 
plotted  in  figure  38.  A  direct  calculation  of  MTF  is  to  take  the  derivation  of  the  ESF  and 
result  a  LSF,  then  take  Fourier  transform.  As  it  can  been  seen,  the  ESF  in  our  case  is 
not  smooth  enough  for  a  direct  calculation.  An  alternative  technique  is  to  use  an  ESF 
fitting  procedure.  In  the  fitting  method,  the  ESF  is  represented  with  a  term  (1-exp)  and 
an  error  function  (erf)  as  in  equation  (10).  MFT  is  calculated  by  equation  (9)  where  a,  b, 
c  and  d  are  fitting  parameters  from  equation  (10).  The  fitting  curve  is  plotted  in  figure  38 
as  the  solid  line. 


ESF(x)  =  a{  1  -  exp {-b  \  x  -  x()  | )}  +  c  •  erf  (dU2  \  x  -  x0 1)  0^) 

yTr{f)  c  ■  exp(-;rz/2  /  d)  +  a(]  +  4k2  f 2  /  b2)~'  (11) 

(c  +  cz) 


The  calculated  MTF  was  compared  with  that  from  scanned  optics  in  figure  40. The  MTF 
from  scanned  optics  was  calculated  from  a  slit  image.  Optic  structure  was  smeared  out 
with  scanning,  so  no  background  deduction  is  necessary.  An  ideal  MTF  was  also 
plotted  in  figure  40.  An  ideal  MTF  is  the  one  when  we  assume  the  lens  introduces  only 
the  1.86  factor  magnification  and  no  image  degradation.  It  is  measured  with  a  slit  image 
without  optics,  and  its  result  was  multiplied  by  the  magnification  of  the  capillary  lens. 
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The  measured  MTF  with  lens  should  no  better  than  the  ideal  MTF.  However,  for  some 
reason,  the  MTF  for  stationary  lens  is  better  than  the  ideal  MTF.  It  is  possible  that  the 
background  elimination  process  failed  the  MTF  measurement.  Another  possibility  is  that 
the  magnification  of  the  lens  is  different  from  1.86,  since  the  magnification  was 
measured  separately  at  a  different  time. 

5.3.4.  Simulations 

Modeling  defects,  including  waviness,  bending,  and  blockages,  require  sophisticated 
simulations.  The  simulations  plotted  in  Figure  25  and  Figure  26  used  the  same 
roughness  and  waviness  correction  as  that  used  for  fiber  D  in  Figure  4.  A  linear 
tapered  lens  is  ideal  for  mammography,  because,  as  shown  above,  any  curvature  can 
cause  transmission  reductions.  However,  the  lenses  currently  available  have  nonlinear 
profiles.  The  profile  is  not  circular,  but  to  simplify  the  simulation,  we  assume  that  the 
profile  is  circular.  A  uniform  bending  will  provide  the  best  transmission  performance  for 
a  given  bending  angle,  so  the  simulation  will  overestimate  the  lens  performance.  A 
further  simplification  assumes  the  lens  has  a  fixed  channel  size,  although  the  channels 
actually  taper  with  the  lens. 

The  geometry  for  uniformly  bent  nonlinear  lenses  is  shown  in  figure  44. The 
simulation  method  for  a  uniformly  bent  lens  is  an  extension  of  the  simulation  for 
uniformly  bent  polycapillary  fibers.  The  simulation  is  done  by  sampling  channels  in  the 
input  cross  section  of  the  lens  with  an  even  step  size.  In  the  simulation  for  a  bent  fiber  it 
is  assumed  that  the  capillary  bends  in  -y  direction,  and  z  direction  is  along  the  input  axis 
of  the  fiber  In  order  to  reuse  the  old  simulation  code,  the  coordinate  of  the  simulation  x- 
y-z  is  rotated  to  x’-y’-z’,  where  z’  is  the  direction  of  channel  axis  at  the  input  end  and  -y’ 
is  the  channel  bending  direction.  Assuming  the  source  position  and  the  capillary 
position  at  input  end  is  (xs\  ys’,  zs’)  and  (0,  0,  z\)  respectively,  then  the  source  capillary 
distance  is  z(-  zs’.  The  origin  of  the  coordinate  is  selected  to  be  the  focus  point  of  the 
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Figure42.  Transmission  of final  half  of  A  (A2),  and 
its  two  constituent  pieces,  the  first  5.5  cm,  A21  and 
the  final  5.5  cm,  A22.  A21  contains  a  defect. 


Energy  (keV) 

Figure43.  Transmission  of  original  optic  compared 
to  two  pieces.  B  is  the  first  10  cm  (closest  to  the 
source),  A  is  the  final  22  cm. 
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lens. 

Another  parameter 
needed  by  the  simulation  is  the 
bending  curvature  of  the 
selected  channel.  This  value 
varies  with  the  position  of  the 
channel.  The  calculation  of  the 
bending  curvature  is 
demonstrated  in  figure  45. 

First  it  is  assumed  that  the 
channels  are  evenly  distributed 
on  both  input  and  output  cross 
section.  Let  input  and  output  Figure45.  Curvature  estimation  for  lens  simulation. 
diameter  of  the  lens  be  n  and  r0, 

and  focal  distance  oH  be  f ,  and  the  lens  length  He  be  L,  then 

~AB  =  ~Jxf  +  y~2  ■  CE  =  TB*rJr: 
f  t  (f  +  L)=lB  tCD •  d  =  CE-CD 

For  a  small  bending  angle, 

a*L/  R  =  d/ (LI  2) 

where  R  is  the  bending  radius,  d  is  the  distance  between  point  D  and  E  as  shown  in 
Figure,  and  a  is  the  bending  angle.  Thus 

R  =  L2  /  2d  (12) 


Figure  44.  Geometry  for  lens  simulation 

The  simulated  results  in  Figure  26  and  Figure  27  show  higher  transmission  than 
the  experimental  results.  This  indicates  that  the  lens  has  more  reduction  in  transmission 
due  to  the  bending.  This  is  expected,  because  the  lens  does  not  have  a  uniform 
bending. 
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5.3.5  Potential  for  Lead  Glass 

Transmission  from  early  fibers  was  poor  in  the  whole  energy  range.  Recent 
fibers  are  more  promising.  Transmission  for  an  early  95  mm  long  and  a  more  recent 
125  mm  long  lead  glass  polycapillary  fiber  is  shown  in  Figure.  The  transmission  drop  at 
lower  energies  may  be  caused  by  partially  blocked  channels.  Transmission  can  also  be 
improved  by  using  a  short  optic.  Almost  total  scatter  rejection  is  expected  for  a  lead 
glass  optic  longer  than  30  mm. 


5.3.5. 1  Calculation  of  Contrast  Improvement  and  Relative  SNR 
Lenses  with  large  diameters  are  not  currently  available.  Their  potential  behavior  can  be 
predicted  by  scaling  up  existing  measurements  to  a  large  area  lens,  as  shown  in  figure 
47.  The  value  of  h  is  the  length  of  the  optic.  For  a  point  P  on  the  film,  scattered 
radiation  could  come  from  any  direction  in  the  solid  angle  Qpatient>  which  is  the  solid 
angle  subtended  by  the  patient.  The  transmission  for  scattered  x  rays  can  be  calculated 
from  the  mass  aborption  coefficients.  The  calculated  scatter  transmission  depends  on 
the  length  of  the  photon  path  through  the  optic.  Making  a  very  conservative  estimation, 
we  can  fix  this  length  to  the  length  of  the  optic.  There  is  also  a  small  part  of  the 
scattered  x  rays  that  have  almost  the  same  direction  as  those  primary  x  rays.  This  part 
of  the  scattered  x  rays  with  an  incident  angle  less  than  critical  angle  could  still  pass 
through  the  channel.  So  the  total  scatter  transmission  of  the  lens  is 


T 


N.  Tn  ■  nd 

S— optic  rjr>  P  c 


s-optic 


N. 


Q 


patient 


(13) 


where  Ns.0ptic  is  the  number  of  scattered  x  rays  with  the  optic,  Ns  is  the  number  of 
scattered  x  rays  without  the  optic,  Tsis  the  transmission  for  large  angle  scattered  x  rays, 
Tp  is  the  transmission  for  primary  x  rays,  Qpatient  is  the  angle  subtended  by  the  patient,  0C 
is  the  critical  angle  for  transmission  in  the  channels.  Because  the  critical  angle  is  on  the 
order  of  10 3  radians,  the  second  term  is  very  small,  and  can  be  ignored.  So  scatter 
transmission  of  the  lens  can  be  estimated  to  be  equal  to  large  angle  scatter 
transmission,  Ts,  which  can  be  calculated  from  equation  19  and  nearly  scatter  free 
lenses  can  be  obtained  by  using  a  suitable  length  optic.  In  the  above  equation,  we 
assumed  that  the  angular  distribution  of  the 
scattered  x  rays  is  isotropic.  If  the  scatter 
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Figure  46.  Transmission  as  a  function  of 

photon  energy  for  a  lead  glass  capillary.  Figure  47.  An  ideal  large  area  lens. 

The  fibers  are  105-120  mm  in  length. 
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Table  IV.  Transmissions  for  three  lenses.  The  results  for  the  distribution  is 

not  isotropic, 
the  second 
term  in 
equation  (13) 
can  be  larger 
than  expected, 
but  will  still  be 
very  small,  and 
this  assumption 
will  not 
significantly 
effect  the 
result.  In 
equation  (13) 

Ts  ,  Tp  and  6C 
are  all  energy 
dependent.  In 
the  case  where 
x  rays  with  a 
wide  x-ray 
spectrum  are 
used,  Ts  and  Tp 

should  be  replaced  by  the  average  large  angle  transmission  and  primary  transmission. 
To  avoid  underestimating  scatter  transmission,  6C  can  be  set  to  be  the  critical  angle  for 
the  lowest  energy  in  the  spectrum. 

The  scatter  and  primary  transmission  of  the  optic  determine  the  performance  of  the 
capillary  lens.  For  example,  if  the  scatter  fraction  before  the  capillary  lens  is  F,  then  the 
scatter  fraction  after  the  lens,  Foptic,  is 


ideal  borosilicate  lens  are  measured  results  for  the  center  part  of 
the  prototype  lens.  The  results  for  lead  glass  lens  are  simulated 
with  120  kVp  spectrum  and  a  W  target. _ 


Lens 

# 

Lens  type 

Length 

(mm) 

Photon 

Energy 

Primary 

transmission 

T  p-ODtiC 

Scatter 

transmission 

T s-optic 

1 

Ideal 

borosilicate 
Lens.  Scaled 
up  from  the 
center  part  of 
the  prototype 
lens 

166 

20  KeV 

0.6 

0.0 

45  KeV 

0.55 

0.002 

60  KeV 

0.55 

0.005 

70  KeV 

0.45 

0.03 

2 

Actual 
Borosilicate 
prototype  lens 

166 

20  KeV 

0.33 

0.0 

3 

Lead  glass 

30 

69  KeV 

0.63 

0.0076 

F  ,  = - 

Optic  p  '  p 


F-T 

s-optic 


KF 


s -optic 


+  (1  -F)-T, 


p -optic 


1  -  F(1  -  K) 


(14) 


where  Tp.0ptic  is  the  primary  transmission  of  the  lens  and  K  is  the  ratio 

T„ 

K  = 


x— optic 


(15) 


p-optic 


The  contrast  enhancement  achieved  by  using  the  optics  is 

Coplic  =  \-F„plic  =  1  (16) 

Cno.op,,  1  -F  FK  +  (1  -  F)  ' 


The  contrast  enhancement  alone  is  not  the  only  relevant  quality  factor.  The  signal  to 
noise  ratio  (SNR)  is  another  important  quality  factor.  According  to  the  Rose  model, 

k2  =  C2(t>A  ,  (17)  (17) 

where  k  is  the  SNR,  C  is  the  contrast,  ®  is  the  photon  flux,  and  A  is  the  area  of  the 
target.  In  order  to  compare  the  SNR  for  different  scatter  rejection  devices,  relative  SNR 
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Scatter  fraction 


was  used.  Relative  SNR  is  the  SNR 
normalized  by  the  SNR  without  any 
scatter  rejection  device, 

kr=^j  =  ~{FTs+(\~F)Tp),  (18) 

*0  '-'0 

where  kr  is  the  relative  SNR.  C  and 
Co  are  the  contrast  with  and  without 
out  the  scatter  rejection  device, 
respectively.  Ts  is  the  scatter 
transmission  and  Tp  is  the  primary 
transmission  of  the  device. 


Figure  48.  Contrast  improvement  vs.  scatter  fraction  for 
the  three  lenses.  The  results  were  compared  with  the 
performance  of  two  commercial  grids. 

5.3. 5.2  Performance  of  T hree  Hypothetical  Lenses 

Performance  estimates  of  three  hypothetical  lenses  are  given  in  Table  IV.  The  ideal 
borosilicate  lens  is  scaled  up  from  the  center  part  of  the  measured  lens.  Simulated 
transmission  was  used  for  the  lead  glass  lens.  The  performance  of  the  three  lenses 
was  compared  with  the  performance  of  the  conventional  grids  described  in  Table  V. 

The  primary  and  scattered  transmissions  of  the  different  grids  in  Table  V  were 
measured  at  69  keV,  which  is  the  average  energy  of  the  120  kVp  tungsten  spectrum, 
and  with  scatter  fraction  of  85%,  by  C.  E.  Dick  and  J.  W.  Motz.  It  is  assumed  that  the 
scatter  and  primary  transmission  are  the  same  under  different  scatter  fraction.  This 
assumption  is  correct  if  the  angular  distribution  of  the  scattered  x  rays  does  not  change. 
The  performance  of  the  grids  is  assumed  to  be  energy  independent.  They  might 
perform  better  at  lower  photon  energies  if  the  scatter  transmission  is  due  to  the  lead 
plates  which  are  not  thick  enough  to  absorb  the  incident  x  rays.  However,  it  was 
assumed  the  dominant  mechanism  for  scatter  transmission  was  x  rays  passing  through 
the  interstitial  material.  The  contrast  enhancement  and  relative  SNR  were  calculated 

under  different  scatter  fractions  and  plotted  in 

Table  V.  Performance  for  several  figures  48  and  49. 


commercial  grids,  meassured  at  a 
scatter  fraction  of  85%. 18 


Grid  ratio 

Primary 
transmission 
at  69  keV 

Scatter 
transmission 
at  69  keV 

6 

0.73 

0.27 

8 

0.64 

0.13 

12 

0.61 

0.071 
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For  most  medical  x-ray  imaging  applications,  the  effective  energy  ranges  from  20  to  70 
keV.  Mammography  has  an  effective  energy  of  20  keV.  The  results  in  figure  48  show 
that  contrast  improvement  of  all  the  three  hypothetical  lenses  are  better  than  or  as  good 
as  that  of  the  three  grids.  Lens  #2  is  scaled  up  from  the  whole  measured  lens.  This  lens 
has  good  contrast  improvement  because  of  its  total  scatter  rejection.  However  its 
relative  SNR  is  lower  than  that  of  the  grids  when  the  scatter  fraction  is  less  than  80% 
because  of  its  low  primary  transmission.  The  primary  transmission  of  this  lens  is  only 
33%  and  the  primary  transmissions  for  all  the  three  commercial  grids  are  better  than 
60%.  So  with  this  lens  more  exposure  will  be  necessary  to  the  patient.  Overall  primary 
transmission  has  to  be  improved  for  practical  use.  Lens#1  is  scaled  up  from  the  center 
part  of  the  measured  lens.  From  the  calculations,  summarized  in  figures  48  and  49,  its 
performance  is  significantly  better  than  those  of  the  commercial  grids  for  effective 
energies  of  20-45  keV.  This  lens  will  not  only  improve  the  contrast  in  the  imaging  but 
also  reduce  the  exposure  time  to  the  patient.  So  it  is  ideal  for  mammography. 

Assuming  that  the  lead  glass  capillaries  have  the  same  quality  as  that  of 
borosilicate  fibers,  we  can  estimate  the  performance  of  a  lead  glass  lens.  Simulations  of 
transmission  spectra  were  performed  for  a  source  to  lens  distance  of  1 .5  m,  and  the 
source  operated  at  120  kVp  with  a  0.5  mm  spot  size.  For  a  perfect  lens,  each  channel 
in  the  lens  is  tapered.  With  a  tapered  channel,  fewer  average  reflections  will  be  needed 
for  the  photons  to  pass.  So  theoretically,  transmission  with  a  perfect  lens  is  larger  than 
for  a  straight  capillary  if  they  have  the  same  quality.  The  transmission  is  quite  sensitive 
to  the  capillary  length.  Therefore,  as  long  as  the  lens  has  low  enough  scatter 


transmission,  the  length  of  the 
lens  should  be  kept  as  short  as 
possible. 

The  performance  estimation 
results  of  the  lead  glass  lens  in 
figure  49  are  based  on  the 
simulated  transmission  with 
length  of  30  mm.  The  scatter 
transmission  of  this  lens  is  less 
than  1%. 
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Scatter  fraction 


Figure  49.  Relative  Signal  to  Noise(SNR)  vs.  scatter  fraction  for 
three  hypothetical  lenses.  The  results  were  compared  with  the 
performance  of  three  commercial  grids. 


Primary  transmission  as  high  as  that  of  grid  is  also  obtained  with  simulation.  Compared 
to  the  grid  performance,  this  lens  has  higher  contrast  enhancement  and  relative  SNR. 
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So  this  lens  will  not  only  improve  the  contrast  of  the  image,  but  also  reduce  the 
radiation  dose  to  the  patient. 


5.3.5.3.  Leaded  Glass  Single  Fiber  Measurements  and  Analysis 

The  experimental  apparatus  for  leaded  glass  single  fibers  measurements  is  shown 
in  figures  2  and  3.  Description  of  five  types  fibers  measured  and  their  best  fit  simulation 
parameters  are  listed  in  Table  VI. 

Figure  50  shows  a  source  scan  plot  of  transmission  of  fiber  B60  as  a  function  of 
source  location  for  eight  different  energy  windows  of  approximately  2  keV  in  width.  The 
transmission  plots  are  narrower  for  the  higher  energy  windows,  as  expected  due  to  the 
decrease  in  critical  angle  with  photon  energy.  For  leaded  glass  fibers,  because  of  the 
high  absorption  of  lead,  the  background  signal  from  photons  with  incidence  angle 
greater  than  the  critical  angle  is  very  small.  Most  of  these  photons  are  absorbed  rather 
than  cutting  through  the  glass  wall.  The  cut  through  transmission  is  3.6%  at  80  keV  for 
a  30-mm-long  type  B  fiber  and  0.18%  for  a  60-mm-long  fiber.  Measured  and  theoretical 
high-angle  transmission  at  several  energies  are  shown  in  figure  51.  The  theoretical 
values  are  calculated  as 

Tc  =  exp  {-L  ( 1-f )  pp }  (19) 


where  L  is  the  fiber  length,  f  is  the  fractional  open  area  of  the  fiber,  p  is  the  density  of 
the  glass,  and  p  is  the  mass  absorption  coefficient  of  the  glass  obtained  from  the 
tabulated  values.19  The  transmission  results  which  follow  are  channel  transmissions 
from  which  the  very  small  high  angle  cut  through  and  fluorescence  background  have 
been  subtracted. 

Transmission  for  leaded  glass  fibers  listed  in  Table  VI  was  measured  in  the  energy 
range  9-80  keV.  The  measured  transmission  as  a  function  of  energy  along  with  the 


9  keV 


Source  Position  (mm) 


Figure  50.  Normalized  intensity  vs  source 
position,  taken  by  moving  the  source  in  the 
horizontal  direction  perpendicular  to  the  fiber 
with  fiber  and  pinhole  fixed. 


Figure  51.  Large  angle  (cut  through)  transmission 
versus  energy,  measured  with  the  source  moved 
away  from  the  aligned  position  to  an  angle  larger 
than  the  critical  angle.  Dash  lines  are  theoretical 
calculation  using  equation  19. 
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simulation  of  different  bending  radii  for  type  B60  fiber  is  shown  in  figure  52.  A  slight 
bending  can  dramatically  reduce  the  transmission  at  high  energies.  The  simulation  with 
bending  alone  can  not  fit  the  experimental  data  since  at  smallest  radius  the 
transmission  is  too  high  at  lower  energies  and  too  low  at  high  energies.  Bending  has 
very  small  effect  on  the  transmission  at  low  energies.  Since  any  waviness  correction  will 
reduce  the  transmission,  the  radius  of  bend  has  to  be  larger  than  32  m. 


Table  VI.  Fiber  description  and  simulation  parameters 


Fiber 

Simulation 

Fiber 

type 

Outer 

diameter 

mm 

Channel 

size 

pm 

Open 

area 

% 

Length 

mm 

Bending 

R 

m 

Waviness 

w 

mrad 

Layer 

t 

pm 

Roughness 

z 

nm 

■ 

12 

60 

60 

40 

0.26 

45 

1.8 

12 

60 

30 

30 

0.15 

18 

1.7 

12 

0.15 

33 

1.7 

11 

27 

0.19 

33 

1.8 

D 

0.53 

11 

50 

60 

28 

0.20 

43 

1.8 

E 

0.53 

11 

40 

60 

45 

0.28 

47 

2.1 

Figure  52.  Simulated  transmission  spectra  of  a 
type  B60  fiber  for  different  bending  radii  are 
compared  with  the  experimental  data.  The  figure 
shows  that  bending  alone  can  not  fit  the  data 
since  the  smallest  bend  underestimates  the 
transmission  at  high  energies. 
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Energy  (keV) 

Figure  53.  Simulated  transmission  spectra  are 
compared  with  the  experimental  data.  The 
waviness  correction  causes  the  transmission  to 
drop  primarily  in  the  middle  energy  range. 
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Simulations  with  different  waviness  corrections  are  compared  to  the  experimental 
data  of  transmission  spectra  in  figure  53.  The  waviness  alone  can  not  fit  the  data.  The 
simulation  shows  that  the  waviness  causes  the  transmission  to  drop  rapidly  at  medium 
energies  but  not  rapidly  enough  at  higher  energies.  Simulated  transmission  source 
scans  with  different  waviness  are  compared  to  the  experimental  data  at  9  keV  in  figure 
54.  The  simulation  with  a  waviness  of  0.5  mrad  fits  the  spectrum  source  scan  fairly  well 
at  9  keV.  However,  that  was  seen  to  be  too  large  for  the  transmission  spectrum  of  figure 
53.  Figure  55  compares  the  experimental  data  to  simulations  at  40  keV,  and  shows  that 
the  maximum  value  of  waviness  is  0.15  mrad. 


9  keV :  ■  exp 


Source  Position(mm) 
Figure  54.  Simulated  source  scan  curves 
with  different  waviness  are  compared  with 
the  experimental  data  at  9  keV.  A  waviness 
of  0.5  mrad,  which  is  too  large  at  higher 
energies,  is  needed  to  fit  the  data. 


40  keV  ■  exp 


Figure  55.  Simulated  source  scan  curves  with 
different  waviness  are  compared  with  the 
experimental  data  at  40  keV.  Since  the  simulated 
curve  with  a  waviness  of  0.15  mrad  fits  the  data 
fairly  well,  it  is  used  as  the  value  for  waviness  in 
the  best  fit  simulation  parameters  in  Table  VI. 


9  keV :  ■  exp 


Figure  56.  Simulated  source  scan  curves 
for  a  type  B60  fiber  for  different  roughness 
are  compared  with  experimental  data  at  9 
keV.  The  curve  is  too  low  at  large 
displacements  for  the  highest  value  of 
roughness. 


The  roughness  has  only  a  fairly  small  effect  on 
the  transmission  spectrum.  The  simulated 
transmission  source  scan  curves  with  different 
roughnesses  compared  with  the  experimental 
data  at  9  keV  is  shown  in  figure  56.  The 
transmission  at  high  angles  is  sensitive  to 
roughness.  The  transmission  is  too  small  if  the 
roughness  is  larger  than  1.7  nm. 

By  modeling  channel  blockage  with  a  glass 
filter  layer  of  thickness  t,  the  simulated 
transmission  spectrum  with  the  glass  layer  alone 
fits  the  data  at  low  energies.  The  simulated 
transmission  spectra  with  different  layer 
thicknesses  are  compared  with  the  experimental 
data  in  figure  57.  The  simulated  transmission 
with  best  fit  simulation  parameters,  which  with 
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and  without  including  glass  inclusions  with  a  thickness  of  33  pm  compared  with  the  data 
for  type  B60  fiber  is  shown  in  figure  58.  As  shown  in  figure  58,  the  inclusion  of  a  33  pm 
thick  leaded  glass  layer  gives  a  good  fit  with  the  experimental  transmission  spectrum. 

The  best  fit  simulation  compared  to  the  experimental  source  scan  curves  are  shown 
in  figures  59-62.  The  simulated  transmission  source  scan  curves  fit  the  data  fairly  well 
at  energies  from  20  to  80  keV.  The  source  scan  at  9keV  does  not  fit  the  simulated 
transmission,  as  shown  in  figure  59.  This  may  be  partly  because  the  actual  value  of  the 
transmission  at  9  keV,  as  shown  in  figure  58,  is  quite  low. 


Energy  (keV) 


Figure  57.  Simulated  transmission  with 
different  layer  thicknesses  alone  compared 
with  experimental  data.  The  figure  shows  that 
the  maximum  layer  thickness  for  a  type  B60 
fiber  is  35  pm  to  fit  the  low  energy 
transmission. 
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Figure  58.  Simulated  transmissions  with  and 
without  leaded  glass  filter  for  type  B60  fiber 
are  compared  with  the  experimental  data. 
Best  fit  is  simulated  transmission  with  R=56 
m,  w=0.15  mrad,  z-1. 7  nm,  and  t=33  pm. 


The  primary  effect  of  bending  is  to  reduce  the  transmission  at  high  energies.  Channel 
blockage  reduces  the  transmission  at  low  energies.  Roughness  reduces  the  width  of 
source  scan  at  low  energies.  Waviness  primarily  reduces  the  transmission  and  the 
source  scan  width  at  mid  range  energies. 


Figure  59.  Simulation  of  source  scan  curves 
with  z=1. 7  nm  ,  R=56  m,  and  w=0. 15  mrad, 
and  with  roughness  z=1.7  nm  alone, 
compared  with  the  experimental  data.  Since 
the  actual  transmission  at  this  energy  is 
quite  low,  the  source  scan  curve  might  be 
somewhat  distorted,  and  thus  the 
simulation  does  not  fit  the  data  well. 


•  exp  20  keV 


Source  Position 

Figure  60.  Simulated  scan  curves  compared 
with  the  experimental  data  at  three  photon 
energies. 
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50  k  eV:  ♦  exp 


Figure  61.  The  simulation  of  scan  curves  with 
best  fit  parameters  at  50  and  70  keV, 
compared  with  the  experimental  data. 


C 


Figure  62.  The  simulation  of  scan  curves  with 
best  fit  parameters  at  60  and  80  keV, 
compared  with  the  experimental  data. 


The  best  fit  simulated  transmissions  along  with  the  experimental  data  are  shown  in 
figures  50  and  51 . 20,21  The  outer  diameter  of  the  fibers  are  in  the  0.51-  0.53  mm  range. 
The  fractional  open  area  of  the  fibers  range  from  60%  to  40%.These  fibers  are  thin  and 
flexible,  and  therefore  are  difficult  to  keep  straight  in  the  measurement  apparatus, 
resulting  in  sharp  bending  radii  and  relatively  poor  high-energy  transmission. 

Fiber  type  B  has  high  open  area,  low  channel  blockage,  and  the  least  bending  and 
waviness.  Thus  this  fiber  has  better  transmission  than  the  other  fibers.  Figure  63  shows 
that  the  transmission  of  this  fiber  type  with  a  shorter  length  (B30)  is  higher  at  lower 
energies  because  of  its  thinner  blockage  layer.  This  is  consistent  with  a  statistically 
random  model  of  glass  inclusion.  Because  of  its  high  flexibility,  the  shorter  fiber 
experiences  more  bending  than  the  longer  ones,  and  thus  the  transmission  falls  off 
faster  at  higher  energies.  The  shorter  fiber  has  a  higher  scatter  cut  through,  as  shown  in 


A:  ■  exp 


Energy  (keV) 

Figure  63.  Simulation  of  transmission  spectra 
of  fiber  A,  fiber  B60  (B  with  length  of  60  mm) 
and  fiber  B30  (B  with  length  of  30  mm)  with 
their  best  fit  parameters  listed  in  Table  VI, 
compared  with  the  experimental  data. 


Figure  64.  Simulated  Transmission 
spectra  of  fiber  C,  fiber  D  and  fiber  E  with 
their  best  parameters  listed  in  Table  VI, 
compared  with  the  experimental  data. 
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figure  51.  The  different  open  area,  bending,  waviness  and  blocked  channels  cause  the 
differences  in  the  transmission  among  the  fibers.  Fiber  type  C  and  D  have  the  same 
open  area,  but  fiber  type  C  has  thinner  glass  layer,  as  listed  in  Table  VI.  As  a  result, 
fiber  type  C  has  higher  transmission  at  lower  energies  than  fiber  type  D,  as  shown  in 
figure  64. 

The  high  measured  transmission  (>42%)  for  fiber  type  B  shows  that  the  performance 
of  this  leaded  glass  fiber  is  quite  good  in  the  20-50  keV  energy  range.  The  geometric 
simulation,  with  roughness,  waviness,  bending  and  channel  blockage  corrections, 
performs  well  in  explaining  the  measured  transmission  spectra  above  15  keV.  , 
According  to  this  simulation,  bending  is  the  most  harmful  at  higher  energies,  and  glass 
inclusions  that  block  the  channels  are  harmful  at  lower  energies.  For  type  B  fibers,  even 
though  the  roughness  is  high,  the  waviness  is  low  and  the  bending  radius  is  the  largest 
among  the  fibers.  The  high-energy  transmission  performances  could  be  improved  by 
making  the  fiber  diameter  larger  and  therefore  more  rigid  to  minimize  unintentional 
bending.  Low  energy  performance  is  improved  by  using  shorter  fibers.  Because  of  the 
good  absorption  of  lead,  the  high-energy  photons,  which  are  not  reflected,  are 
absorbed  completely  by  60  mm-long  fibers,  even  at  80  keV.  The  high-angle 
transmission  at  80  keV  is  only  0.2%.  A  shorter  fiber,  which  is  30  mm  in  length  has  a 
transmission  efficiency  more  than  40%  at  20  keV  and  scatter  transmission  of  3.6%  at 
80  keV.  This  high-angle  transmission  is  very  much  less  than  the  16%  for  a  longer  (140 
mm)  borosilicate  fiber.  Thus  short  leaded  glass  polycapillary  x-ray  optics  are  very 
promising  for  scatter  rejection  in  mammography  applications. 


6  Key  Research  Accomplishments 

•  Development  of  extensive  modeling  capability  for  polycapillary  optics  that  includes 
realistic  models  for  the  effects  of  profile  defects  and  shows  extremely  good 
agreement  with  measured  data  and  therefore  the  capability  to  predict  performance 
in  new  geometries. 

•  Measured  polycapillary  fiber  transmission  in  excess  of  65%  at  20  keV  for 
borosilicate  glass.  Measured  high  angle  transmission  (which  corresponds  to  scatter 
transmission)  is  lower  than  1%  for  12  cm  long  borosilicate  glass  optics. 

•  Modeling,  measurement  and  defect  studies  have  led  to  better  manufacturing 
capabilities  for  magnifying  tapered  monolithic  optics.  These  optics  have 
demonstrated  the  ability  to  provide  nearly  ideal  scatter  rejection  and  contrast 
enhancement,  while  at  the  same  time  increasing  the  system  MTF  performance  at  all 
spatial  frequencies.  A  multitaper  optic  has  been  produced  and  tested. 

•  Image  analysis  studies  have  provided  the  means  for  image  artifact  analysis  and 
reduction. 

•  Measured  leaded  glass  polycapillary  fiber  transmission  is  in  excess  of  45%  at 

20  keV.  Measured  high  angle  transmission  is  less  than  4%  for  a  30  mm  long  fiber  at 
80  keV. 


Simulations,  Measurements  and  Image  Processing  for  Capillary  Optical  Digital  Mammography 


35 


7  Reportable  Outcomes 
7.1  Manuscripts 

Lei  Wang,  C.  A.  MacDonald,  and  W.  W.  Peppier,  “Performance  of  Polycapillary  Optics 
for  Hard  X-ray  Imaging,”  accepted,  Medical  Physics. 

Lei  Wang,  W.M.  Gibson,  C.A.  MacDonald,  “Potential  of  Polycapillary  X-Ray  Optics  in 
Medical  Imaging  Applications,”  in  C.A.  MacDonald,  K.A.  Goldberg,  J.R.  Maldonado,  A.J. 
Marker  III,  S.P.  Vernon,  eds.,  EUV,  X-ray,  and  Neutron  Optics  and  Sources,  SPIE 
vol.  3767,  1999. 

S.D.  Padiyar,  M.V.  Gubarev,  Hui  Wang,  W.M.  Gibson,  C.A.  MacDonald 
“Characterization  of  Polycapillary  X-Ray  Collimating  Optics”,  in  C.A.  MacDonald,  K.A. 
Goldberg,  J.R.  Maldonado,  A.J.  Marker  III,  S.P.  Vernon,  eds.,  EUV,  X-ray,  and 
Neutron  Optics  and  Sources,  SPIE  vol.  3767,  1999. 

S.D.  Padiyar,  H.  Wang,  W.M.  Gibson,  C.A.  MacDonald,  M.V.  Gubarev,  “Beam 
Collimation  Using  Polycapillary  X-Ray  Optics  For  Large  Area  Diffraction  Applications,” 
in  Advances  in  X-ray  Analysis,  43,  Proceedings  of  the  48th  Denver  X-ray  Conference, 
1999. 

Hui  Wang,  Lei  Wang,  W.M.  Gibson,  C.A.  MacDonald,  “Simulation  Study  of  Polycapillary 
X-Ray  Optics,”  in  X-Ray  Optics,  Instruments,  and  Missions,  R.B.  Hoover  and  A  B  C. 
Walker  II,  eds.  SPIE  Vol  3444,  pp.  643-651,  July  1998. 

B. K.  Rath,  W.M.  Gibson,  Lei  Wang,  B.E.  Homan  and  C.A.  MacDonald,  “Measurement 
and  Analysis  of  Radiation  Effects  in  Polycapillary  X-ray  Optics,”  Journal  of  Applied 
Physics,  83,  no. 12,  pp.  7424-7435,  June  15  1998. 

Suparmi,Cari,  W.M.  Gibson,  C.A.  MacDonald,  "Development  of  Polycapillary  X-Ray 
Optics  for  Scatter  Rejection,”  in  Advances  in  Laboratory-Based  X-Ray  Sources  and 
Optics,  SPIE  Vol.  4144,  2000. 

Suparmi,  Cari,  Lei  Wang,  Hui  Wang,  W.M.  Gibson,  C.A.  MacDonald,  ’’Measurement 
and  Analysis  of  Leaded  Glass  Polycapillary  Optic  Performance  for  Hard  X  Rays,”  to  be 
submitted  for  Journal  of  Applied  Physics. 

Cari,  Suparmi,  S.D.  Padiyar,  W.M.  Gibson,  C.A.  MacDonald,  C.D.  Alexander,  M.K.  Joy, 

C. H.  Russel,  Z.W.  Chen,  “Characterazation  of  a  Long  Focal  Length  Polycapillary  Optic 
for  High  Energy  X  Rays,”  in  Advances  in  Laboratory-Based  X-Ray  Sources  and 
Optics,  SPIE  Vol.  4144,  2000. 


6imulations,  Measurements  and  Image  Processing  for  Capillary  Optical  Digital  Mammography 


36 


7.2  Presentations 

Lei  Wang,  W.M.  Gibson,  C.A.  MacDonald,  “Potential  of  Polycapillary  X-Ray  Optics  in 
Medical  Imaging  Applications,”  SPIE  1999. 

Suparmi,  Cari,  W.M.  Gibson,  C.A.  MacDonald,  "Development  of  Anty  Scatter  Optics 
Using  Leaded  Glass  Polycapillary, "  SPIE  2000. 

7.3  Poster  Presentations 

Hui  Wang,  Lei  Wang,  W.M.  Gibson,  C.A.  MacDonald,  “Simulation  Study  of  Polycapillary 
X-ray  Optics,”  SPIE  San  Diego  7/98. 

Suparmi,  Cari,  W.M.  Gibson,  C.A.  MacDonald,  "Measurement  and  Analysis  of  Leaded 
Glass  Polycapillary  Optics  for  Mammography,"  Era  Hope  Meetings,  Atlanta,  2000. 

7.4  Employment/Research  Received 
7.4.1  Research  Employment 

Lei  Wang,  Postdoctoral  Residency,  University  of  California  at  San  Francisco,  1999. 


8  Conclusions 

Some  manufacturing  problems  have  occurred  with  scaling  up  the  optics  to 
clinical  size.  These  problems  have  been  addressed  by  the  studies  of  defect  properties. 
Notwithstanding  the  problems,  the  optics  tested  have  performed  extremely  well. 
Multifiber  collimating  optics  and  monolithic  linear  magnifying  tapers  have  been  studied. 
Scatter  rejection  is  very  high  from  all  the  optics.  The  linear  tapers  resulted  in  nearly 
ideal  contrast  enhancements  and  simultaneously  increase  in  MTF  at  all  spatial 
frequencies.  This  will  result  in  an  improvement  in  resolution  regardless  the  inherent 
spatial  resolution  of  the  detector.  The  rapidly  growing  modeling  capability  has  lead  to  a 
real  development  of  understanding  of  the  nature  of  polycapillary  defects,  which  has 
lead  to  improvements  in  the  manufacturing  processes.  Further,  the  increased 
confidence  in  the  modeling  allows  future  lens  geometries  to  be  effectively  “tested”  in 
simulation  so  that  design  parameters  can  be  rapidly  optimized.  Short  leaded  glass 
polycapillary  optics  which  have  good  transmission  at  mammographic  energies  provide 
excellent  for  scatter  rejection. 
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1.  ABSTRACT 

The  properties  of  borosilicate  glass  single  polycapillary  optics  have  been  previously 
studied  in  10-80  keV  energy  range.  However,  many  medical  imaging  and  industrial 
applications  for  x  rays  require  large  area  optics  with  good  scatter  rejection.  Since  shorter 
optics  are  easier  to  manufacture,  optics  with  a  shorter  length  would  provide  a  faster  route 
to  bringing  the  benefit  of  polycapillary  optics  to  these  applications.  Leaded  glass  would 
allow  the  optic  to  be  much  shorter  and  still  give  good  contrast  enhancement,  because  of 
the  superior  absorption  of  leaded  glass.  In  order  to  investigate  the  feasibility  of  using 
leaded  glass  polycapillary  optics  for  these  applications,  measurements  and  simulations 
have  been  performed  on  the  behavior  of  leaded  glass  polycapillary  fibers  in  the  9-80  keV 
energy  range.  These  fibers  have  outer  diameter  about  0.5  mm,  channel  diameter  of  12  pm 
and  open  area  from  40%  to  60%.  The  transmission  efficiencies  of  these  fibers  of  different 
types  and  lengths  were  measured  as  a  function  of  source  location  and  x-ray  energy.  The 
measurements  were  analyzed  using  a  geometrical  optics  simulation  program,  which 
included  roughness,  waviness,  bending  effects  and  a  leaded  glass  filter  layer.  Despite  low 
transmission  at  low  energies,  leaded  glass  polycapillary  optics  with  a  length  of  30-60  mm 
seem  promising  for  many  high  energy  (>20  keV)  x-ray  applications.  Type  B  fibers 


60  mm  in  length  have  transmission  efficiency  of  50%  in  the  35-40  keV  energy  range, 
which  is  83%  of  fiber’s  open  area,  and  very  low  scatter  transmission  of  0.1 8%  at  80  keV. 
The  same  fiber  with  the  shorter  length  of  30  mm  has  higher  transmission  at  lower 
energies,  and  more  than  52%  in  the  20-  40  keV  energy  range. 
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2.  INTRODUCTION 


Polycapillary  optics,  consisting  of  arrays  of  tiny  hollow  glass  tubes  with  diameters  on 
the  order  of  few  microns,  is  a  relatively  new  technology  for  controlling  x-ray  beams.  A 
polycapillary  fiber  is  a  capillary  bundle  with  hundreds  or  thousands  of  hollow  channels, 
as  shown  in  figure  1.  X  rays  striking  the  interior  of  the  glass  tubes  at  grazing  incidence 
are  guided  down  the  tubes  by  total  external  reflection.  Arrays  of  curved  or  tapered 
capillaries  can  be  used  to  focus,  collimate  and  filter  x-ray  radiatioa  Polycapillary  optics 
have  potential  usefulness  in  many  applications,  including  x-ray  lithography,  medical 
imaging,  crystallography  and  astronomy.1'9  Properties  of  borosilicate  polycapillary  fibers, 
fabricated  optics,  and  their  applications  have  been  studied  extensively  in  the  8-80  keV 
energy  range.10"15  Leaded  glass  fiber  optics  are  advantageous  for  excellent  scatter 
rejection  compared  to  borosilicate  glass  fibers  because  of  their  superior  absorption. 

In  this  paper,  measurement  results  for  variety  of  leaded  glass  single  fibers  in  the 
energy  range  from  9  to  80  keV  are  presented.  The  results  were  analyzed  using  a 
geometrical  optic  simulation  program.  The  use  of  the  simulation  to  analyze  the 
experimental  data  from  60  mm  long  type  B  fibers  is  demonstrated.  Finally,  the  fitting 
processes  were  used  to  determine  the  best  fit  parameters  of  the  simulation  for  five  other 
leaded  glass  fibers.  The  measurements  show  the  potential  for  applications  at  energies 
from  25  to  80  keV  for  short  length  optics.  The  analyses  from  the  simulation  indicate  that 
waviness  and  bending  are  harmful  to  the  transmission  at  higher  energies,  and  partially 
blocked  channels  are  more  harmful  for  lower  energies. 
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3.  BASIC  PRINCIPLE 


Capillary  x-ray  optics  use  multiple  total  external  reflection  to  guide  grazing  incidence 
x  rays  in  arrays  of  tiny  hollow  glass  tubes.  X-ray  photon  energies  are  much  larger  than 
the  plasma  energies  of  glasses,  which  are  tens  of  electron  volts.  In  this  regime,  the  real 
part  of  refractive  index  of  glass  can  be  simply  approximated  by16 

n2  =s  / £0  =  l-(a?p  / a? )  0) 

where  n  is  the  index  of  refraction  of  glass,  s  is  the  dielectric  constant  of  the  glass,  <x>  is 
the  x-ray  frequency,  (Op  is  the  plasma  frequency  of  the  glass,  and  £<?is  the  dielectric 
constant  for  vacuum.  The  plasma  frequency  (Op  for  leaded  glass  is  small  compared  to  the 
photon  frequency,  so  that  n  is  slightly  less  than  unity.  Therefore,  total  external  reflection 
can  occur  at  the  air  to  glass  interface.  Using  Snell’s  Law,  the  critical  grazing  angle  0C  can 
be  found  from 

sin(~  -0c)  =  nx  sin(  j).  (2) 

Therefore, 

00  =  (Op/  (o  (3) 

For  leaded  glass  capillaries  that  were  used  in  the  experiments  described  here,  the  critical 
angle  for  reflection  is 

0c-35  (mrad)  /E  (keV)  (4) 

The  higher  the  photon  energy,  the  smaller  is  the  critical  angle. 
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4.  GEOMETRIC  SIMULATION 


In  the  geometrical  optic  simulation  used  in  this  article,  x-ray  transmission  through 
hollow  glass  fibers  is  simulated  by  tracing  a  large  number  of  x  rays  through  fibers.  The 
simulation  includes  four  fiber  quality  defects:  roughness;  waviness;  channel  blockage;  and 
also  profile  corrections  that  can  be  approximated  by  uniform  bending.  It  had  been  shown 
previously  that  geometric  simulations  with  these  four  best  parameters  are  in  good 
agreement  with  experimental  data  for  borosilicate  single  fibers  at  1-  80keV. 13,17 

The  geometric  algorithm  is  a  two  dimensional  approximation.  Because  of  the  small 
critical  angle,  the  velocity  of  the  x-ray  photon  along  the  capillary  axis  is  nearly  equal  to  c, 
the  velocity  of  light.  The  trajectory  of  an  x-ray  can  be  reduced  to  the  two  dimensional 
motion  inside  the  capillary  cross  section.13  Deviations  of  the  channel  from  a  straight  path 
are  approximated  by  uniform  bending  described  by  a  bending  radius,  R.  The  bending  of  the 
channel  make  the  apparent  motion  of  the  x  ray  within  the  channel  cross  section  similar  to 
that  of  a  classically  accelerated  particle. 

Surface  roughness  is  formally  parameterized  by  a  correlation  length  and  root-mean 
square  displacement  of  rough  surface,  z.  Because  changes  in  the  correlation  length  change 
the  effective  roughness  height,  and  due  to  the  lack  of  actual  correlation  data,  the  correlation 
length  was  fixed  at  6  pm  for  all  simulations.  Changing  this  value  would  cause  a 
proportional  scaling  all  of  reported  roughness  heights.  The  value  was  chosen  for 
borosilicate  fibers  so  that  the  roughness  heights  agree  with  atomic  force  microscopy  data.13 
The  surface  roughness  decreases  the  apparent  reflectivity  of  the  channel  walls  and  therefore 
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the  transmission  in  capillary  channels.  The  third  capillary  surface  quality  parameter  is 
waviness,  which  occurs  on  a  spatial  range  shorter  than  the  capillary  length  and  longer  than 
roughness.  The  average  effect  of  the  waviness  can  be  considered  as  a  random  tilt  of  glass 
wall.  The  maximum  random  tilt  angle  w  is  an  adjustable  parameter.  In  the  simulation,  these 
tilt  angles  are  normally  distributed  with  a  mean  value  of  zero  and  a  standard  deviation  w. 
The  fourth  polycapillary  defect  is  channel  wall  blockage,  which  is  modeled  by  a  filter  layer 
of  glass  of  thickness  t.  The  layer  represents  glass  inclusions  which  randomly  occur  along 
the  fiber. 

5.  EXPERIMENTAL  APPARATUS  AND  TECHNIQUE 
5.1  Experimental  Apparatus 

The  experimental  apparatus  is  shown  in  figure  2. 

5.1.1.  Source 

The  x-ray  generator  used  in  the  experiment  was  a  low  current  Microfocus  MS50 
with  a  50  pm  spot  size,  tungsten  target,  and  a  maximum  operating  voltage  of  100  kVp. 
The  source  head  was  mounted  on  a  stage  which  could  be  moved  in  the  two  directions 
transverse  to  the  x-ray  beam.  In  order  to  reduce  the  background  of  scattered  x  rays,  this 
movable  source  was  enclosed  in  a  6-mm-thick  lead  box  with  a  100x175  mm  aperture  in 
the  front.  The  source  to  fiber  distance  was  1100  mm.  In  the  measurement  for  photons 
with  energy  higher  than  50  keV,  a  5-mm-thick  aluminum  plate  was  used  as  a  filter 
between  the  source  and  the  fiber  to  remove  the  low  energy  photons,  and  reduce  the  dead 
time  of  detector. 
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5.1.2.  Pinhole 


A  200  diameter  pinhole  through  4  mm  of  lead,  2  mm  tungsten,  and  1  mm  of  tantalum 
(attached  together),  was  placed  5mm  away  from  the  fiber  at  the  entrance  end.  The  pinhole 
was  attached  to  a  6-mm-thick  lead  shield.  The  pinhole  and  lead  shield  together  block 
scattered  x  rays  and  leakage  around  the  outside  of  the  fiber.  The  pinhole  was  smaller  than 
the  fiber  but  still  covers  hundreds  or  more  channels.  Its  small  size  was  chosen  not  only  to 
avoid  the  leakage  around  the  fiber  and  reduce  dead  time  in  the  detector,  but  also  to  keep  the 
source-fiber  distance  as  small  as  possible  while  keeping  the  photon  entrance  angle  at  the 
edges  of  the  pinhole  to  less  than  the  critical  angle  at  energies  up  to  80  keV. 

5.1.3.  Fiber 

The  fiber  was  held  straight  by  a  finely  machined  groove  in  an  aluminum  plate  and 
covered  by  iron  powder  to  prevent  x-ray  leakage  around  the  fiber.  The  aluminum  plate 
was  also  mounted  on  a  stage  which  can  be  translated  in  two  orthogonal  directions 
transverse  to  the  beam.  All  stages  were  mounted  on  rail  carriers  which  could  be  moved 
along  the  beam  directions.  Five  kinds  of  fibers,  fisted  in  Table  I,  were  measured.  The 
outer  cross  sections  of  all  five  fibers  are  roughly  hexagonal. 

5.1.4.  Detector 

The  detector  was  a  high  purity  germanium  detector  with  about  200  eV  resolution  at 
5.9  keV  and  550  eV  at  122  keV.  The  detector  was  placed  behind  the  fiber.  The  distance 
between  the  source  and  detector  was  fixed  to  keep  air  absorption  constant.  Motion  control 
and  data  collection  were  all  controlled  with  a  small  computer. 
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6.  MEASUREMENT  RESULTS  AND  ANALYSIS 


6.1.  High  Angle  Transmission 

The  pinhole  was  located  in  the  center  of  the  fiber  to  avoid  leakage.  Figure  3  shows  a 
source  scan  plot  of  transmission  as  a  function  of  source  location  for  seven  different  energy 
windows  of  approximately  2  keV  in  width.  The  transmission  plots  are  narrower  for  the 
higher  energy  windows,  as  expected  due  to  the  decrease  in  critical  angle  with  photon 
energy.  For  leaded  glass  fibers,  because  of  the  high  absorption  of  lead,  the  background 
signal  from  photons  with  incidence  angle  greater  than  the  critical  angle  is  very  small.  These 
photons  are  absorbed,  rather  than  cutting  through  the  glass  wall.  The  cut  through 
transmission  is  3.6%  at  80  keV  for  a  30-mm-long  type  B  fiber  and  0.18%  for  a  60-mm-long 
fiber.  Measured  and  theoretical  high-angle  transmission  at  several  energies  are  shown  in 
figure  4.  The  theoretical  values  are  calculated  using 

T  — ^ ^(1— f)pp 

Tc=e  ^  (5) 

where  l  is  the  fiber  length,  f  is  the  fractional  open  area  of  the  fiber,  p  is  the  density  of  the 
glass,  and  p  is  the  attenuation  coefficient  of  the  glass  obtained  from  the  tabulated  values.18 
The  transmission  results  which  follow  are  channel  transmissions  from  which  the  very  small 
high  angle  cut  through  and  fluorescence  background  have  been  subtracted. 

6.2.  Simulation  Analysis 

Transmission  for  the  leaded  glass  fibers  in  Table  I  was  measured  in  the  energy  of 
9-  80keV.  To  illustrate  the  effect  of  each  simulation  parameter,  the  application  of  the 
simulation  to  the  experimental  data  for  fiber  type  B60  (a  60mm-long  type  B  fiber)  is 
discussed  in  this  section. 
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6.2.1.  Bending  effects 

A  slight  bending  can  dramatically  reduce  the  transmission  at  the  high  energies 
because  of  the  small  critical  angles.  The  measured  transmission  as  a  function  of  energy 
along  with  the  simulations  with  different  bending  radii  for  fiber  type  B60  is  shown  in 
figure  5.  The  simulation  with  bending  alone  can  not  fit  the  experimental  data  since  at 
smallest  radius  the  transmission  is  too  high  at  lower  energies  and  too  low  at  high 
energies.  Since  any  waviness  correction  will  reduce  the  transmission,  the  radius  of 
bending  has  to  be  larger  than  32  m.  Simulations  with  different  bending  radii  and  the 
experimental  transmission  source  scan  at  9  keV  are  shown  in  figure  6.  Even  with  a  sharp 
bend,  the  simulated  transmission  source  scan  at  9  keV  is  not  narrow  enough.  Thus 
additional  corrections  are  necessary. 

6.1.2.  Waviness  effects 

Simulations  with  different  waviness  corrections  are  compared  to  the  experimental 
data  of  transmission  spectra  for  a  type  B60  fiber  in  figure  7.  The  waviness  alone  can  not 
fit  the  data.  The  simulation  shows  that  the  waviness  causes  the  transmission  to  drop 
rapidly  at  medium  energies  but  not  rapidly  enough  at  higher  energies.  Figure  7  shows  that 
the  waviness  correction  has  to  be  smaller  than  0.3  mrad  to  fit  the  data.  Simulated 
transmission  source  scans  with  different  waviness  are  compared  to  the  experimental  data 
at  9  keV  in  figure  8.  The  simulation  with  a  waviness  of  0.5  mrad  fits  the  spectrum  source 
scan  fairly  well  at  9  keV.  However,  that  was  seen  to  be  too  large  for  the  transmission 
spectrum  of  figure  7.  Figure  9  compares  the  experimental  data  to  the  simulations  at  40 
keV,  and  shows  that  the  maximum  value  of  waviness  is  0.15  mrad. 
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6.2.3.  Roughness  effects 

Simulation  of  the  transmission  of  fiber  type  B60  with  and  without  roughness 
corrections  are  compared  to  the  experimental  data  in  figure  10.  For  a  perfect  capillary,  the 
transmission  spectrum  is  flat,  since  the  distance  between  the  source  and  the  fiber  is  large 
enough  to  ensure  that  every  photon  which  hits  the  glass  wall  has  grazing  angle  smaller  than 
its  critical  angle.  The  transmission  is  almost  flat  at  all  energies,  but  for  high  roughness  drops 
more  rapidly  at  medium  energies  than  at  high  energies.  The  transmission  drop  depends  on 
both  the  surface  reflectivity,  which  is  reduced  by  roughness,  and  the  number  of  reflections 
that  the  photon  undergoes  through  the  channel.  The  transmission  is  proportional  to  Rn  where 
R  is  reflectivity,  and  n  is  the  number  of  bounces.  For  a  straight  capillary,  the  average 
number  of  bounces  of  the  photon  along  the  channel  is  small,  usually  less  than  three.  Thus 
the  roughness  has  only  fairly  small  effect  on  the  transmission  spectrum.  The  transmission  as 
a  function  of  source  location  at  9  keV  is  shown  in  figure  11.  Because  photons  experience 
more  than  21  reflections  for  large  source  displacements,  the  transmission  at  high  angles  is 
sensitive  to  roughness  and  becomes  too  small  if  the  roughness  is  larger  than  1.7  nm. 

6.2.4.  Channel  blockage 

The  drop  in  transmission  at  low  energies  can  be  fit  by  modeling  channel  blockage 
with  a  glass  filter  layer  of  thickness  t.  The  simulated  transmission  spectra  with  different 
layer  thicknesses  is  compared  with  the  experimental  data  in  figure  12.  The  simulated 
transmission  spectrum  with  a  glass  layer  thickness  of  35  pm  alone  fits  the  data  type  B60 
fiber  at  low  energies.  To  fit  the  high  energy  transmission  spectrum  data,  finally,  we 
combine  the  waviness  determined  from  source  scan  curve  at  40  keV,  the  roughness 
determined  from  source  scan  curve  at  9  keV,  and  add  the  bending  radius.  To  fit  low 
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energy  transmission,  then  the  simulation  result  must  include  absorption  by  a  leaded  glass 
layer  determined  from  figure  12,  which  must  be  less  than  35  pm.  Figure  13  shows  that 
the  inclusion  of  a  33  pm  thick  glass  layer  gives  a  good  fit  with  the  experimental 
transmission  spectrum. 

The  simulation  with  the  four  best  fit  simulation  parameters  is  compared  to  the 
experimental  source  scans  in  figures  14-17.  The  simulated  transmission  source  scan  at 
medium  and  high  energies  fits  the  data  fairly  well.  The  source  scan  at  9keV  does  not  fit 
the  simulated  transmission,  as  shown  in  figure  14.  This  may  be  partly  because  the  actual 
value  of  the  transmission  at  9  keV,  as  shown  in  figure  13,  is  quite  low. 

6.2.5  Simulation  summary 

The  primary  effect  of  bending  is  to  reduce  the  transmission  at  highest  energies. 
Channel  blockage  reduces  the  transmission  at  lowest  energies.  Roughness  reduces  the 
width  of  source  scan  width  at  these  low  energies.  Waviness  reduces  the  transmission  and 
the  source  scan  at  mid  range  energies. 

6.3.  Explanation  Experiment  Results 

Using  the  fitting  processes  described  in  section  6.2,  the  best  fit  simulation 
parameters  for  fibers  A,  C,  B60,  D  and  E,  which  have  length  60  mm,  and  also  fiber  B30, 
which  has  length  30  mm,  were  obtained.  The  best  fit  simulation  parameters  are  listed  in 
Table  I.  The  best  fit  simulated  transmissions  along  with  the  experimental  data  are  shown 
in  figures  18  and  19.  The  outer  diameter  of  the  fibers  are  in  the  0.51-  0.53  mm  range. 
These  fibers  are  thin  and  flexible,  and  therefore  are  difficult  to  keep  straight  in  the 
measurement  apparatus,  resulting  in  sharp  bending  radii  and  relatively  poor  high-energy 
transmission. 
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Fiber  type  B  has  high  open  area,  low  channel  blockage,  and  the  least  bending  and 
waviness.  Thus  this  fiber  has  better  transmission  than  other  fibers.  Figure  18  shows  that 
the  transmission  of  this  fiber  type  with  the  shorter  length  (B30)  is  higher  at  lower 
energies  because  of  its  thinner  blockage  layer.  This  is  consistent  with  a  statistically 
random  model  of  glass  inclusion.  Because  of  its  high  flexibility,  the  shorter  fiber 
experiences  more  bending  than  the  longer  fiber,  and  thus  the  transmission  falls  off  faster 
at  higher  energies.  The  shorter  fiber  has  a  higher  scatter  cut  through,  as  shown  in  figure  4. 
The  differences  in  the  transmission  among  the  fibers  are  caused  by  the  different  open 
area,  bending,  waviness  and  blocked  channels.  Fiber  type  C  and  D  have  the  same  open 
area,  but  fiber  type  C  has  thinner  glass  layer.  As  a  result,  fiber  type  C  has  higher 
transmission  at  lower  energies  than  fiber  type  D,  as  shown  in  figure  19. 

7.  CONCLUSIONS 

The  high  measured  transmission  (>45%)  for  fiber  type  B  shows  that  the  performance  of 
this  leaded  glass  fiber  B  is  quite  good  in  the  25-50  keV  energy  range.  The  geometric 
simulation,  with  roughness,  waviness,  bending  and  channel  blockage  corrections,  performs 
well  in  explaining  the  measured  transmission  spectra  above  15  keV.  According  to  this 
simulation,  bending  is  the  most  harmful  at  higher  energies,  and  glass  inclusions  that  block 
the  channels  are  harmful  at  lower  energies.  For  type  B  fibers,  even  though  the  roughness  is 
high,  the  waviness  is  low.  The  bending  radius  is  the  largest  among  the  fibers.  The  high- 
energy  transmission  performances  could  be  improved  by  making  the  fibers  larger  in 
diameter  and  therefore  more  rigid  to  minimize  an  unintentionally  bending.  Low  energy 
performance  is  improved  by  using  shorter  fibers.  Because  of  the  good  absorption  of  lead, 
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the  high-energy  photons  which  are  not  reflected  are  absorbed  nearly  completely  by  60  mm- 
long  fibers,  even  at  80  keV.  The  high  angle  transmission  at  80  keV  is  only  0.2%.  At 
energies  below  50  keV,  the  high  angle  transmission  is  negligible  for  the  30  mm  length. 
Thus  short  leaded  glass  polycapillary  x-ray  optics  are  very  promising  for  scatter  rejection 
applications. 
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Table  I.  Fiber  description  and  simulation  parameters 


Fiber 

Simulation 

Fiber 

type 

Outer 

Diameter 
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Channel 

size 

pm 

Open 

area 
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Length 
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Bending  R 
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Waviness  w 
mrad 

Layer  t 
pm 

Roughness  z 
nm 

A 

0.53 

12 

60 

60 

40 

0.26 

45 

1.8 

B30 

0.51 

12 

60 

30 

30 

0.15 

18 

1.7 

B60 

0.51 

12 

60 

60 

56 

0.15 

35 

1.7 

C 

0.51 

11 

50 

60 

27 

0.19 

33 

1.8 

D 

0.53 

11 

50 

. 

60 

28 

0.2 

43 

1.8 

E 

0.53 

11 

40 

60 

45 

0.28 

47 

2.1 
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Detector 


Figure  1.  Photograph  of  a  sample  cross-section  of 
leaded  glass  polycapillary  fiber  with  channel 
diameter  of  12  pm  and  fractional  open  area  of  60%. 


/  » 

Lead  shielding  with 
200  pm  pinhole 

Figure  2.  Experimental  Setup  for 
fiber  measurement 


9  keV 


Source  Position  (mm) 

Figure  3.  Normalized  transmission  vs  source 
position,  taken  by  moving  the  source  in  the 
horizontal  direction  with  fiber  and  pinhole  fixed 


Energy  (keV) 


Figure  4.  Large  angle  (cut  through)  transmission 
versus  energy,  measured  with  the  source  moved 
away  from  the  aligned  position  to  an  angle  larger 
than  the  critical  angle.  Dashed  lines  are 
theoretical  calculations  using  equation  5. 


Energy  (keV) 

Figure  5.  Simulated  transmission  spectra  for 
type  B60  fiber  for  different  bending  radii  are 
compared  with  the  experimental  data.  The  figure 
shows  that  bending  alone  can  not  fit  the  data 
since  the  smallest  bend  underestimates  the 
transmission  at  high  energies. 
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Energy  (keV) 

Figure  7.  Simulated  transmission  spectra  are 
compared  with  the  experimental  data.  The 
waviness  correction  causes  the  transmission  to 
drop  primarily  in  the  middle  energy  range. 


9  keV:  $  exp 


Figure  6.  Simulation  of  source  scan  curves  with 
different  bending  radii  compared  with  the 
experimental  data  at  9  keV.  The  simulated 
curve  does  not  become  narrow  enough  even  for 
a  bend  radius  of  R=20  m. 
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Figure  8.  Simulated  source  scan  curves  F'gure  9-  Simulated  source  scan  curves  with  different 

with  different  waviness  are  compared  with  waviness  are  compared  with  the  experimental  data  at 
the  experimental  data  at  9  keV.  A  wariness  40  keV-  Since  the  simulated  curve  with  a  wariness  of 

of  0.5  mrad,  which  is  too  large  at  higher  °-15  mrad  fits  the  data  fair,y  weU> !t  is  med  as  the 

energies,  is  needed  to  fit  the  data.  vaIue  for  wariness  in  the  best  fit  simulation 

parameters  in  Table  I. 
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Figure  10.  Simulated  transmission  spectra  for 
type  B60  fiber  for  different  roughness  corrections 
are  compared  with  experimental  data 
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Figure  11.  Simulated  source  scan  curves  for  type 
B60  fiber  for  different  roughness  are  compared 
with  experimental  data  at  9  keV.  The  curve  is  too 
low  at  large  displacements  for  the  highest  value  of 
roughness. 
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Figure  12.  Simulated  transmission  spectra  with 
different  layer  thicknesses  alone  compared  with  the 
experimental  data.The  figure  shows  that  the 
maximum  layer  thickness  for  a  type  B60  fiber  is  35 
pm  to  fit  the  low  energy  transmission. 
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Figure  13.  Simulated  transmissions  with  and 
without  a  leaded  glass  filter  layer  for  fiber  B60  are 
compared  with  the  experimental  data.  Best  fit  is 
achieved  with  R=56  m,  w=0.15  mrad,  z=1.7  nm, 
s=6  urn  and  t=33  urn. 


9  keV:  *  exp 


Figure  14.  Simulation  of  source  scan  curves 
with  *=1.7  nm,  R=56  m,  and  w=0.15  mrad, 
and  with  roughness  z=1.7  nm  alone,  compared 
with  the  experimental  data.  The  best  fit 
simulation  does  not  fit  the  data  quite  well, 
possibly  because  the  actual  transmission  at  9 
keV  is  quite  low. 
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Figure  15.  Simulated  scan  curves  with  best 
fit  simulation  parameters  compared  with 
the  experimental  data  at  three  photon 
energies. 
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Figure  16.  The  simulation  of  scan  curves  with 
best  fit  parameters  at  50  and  70  keVs,  compared 
with  th<*  exnpri  mental  data. 
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Figure  17.  The  simulation  of  scan  curves  with  best 
fit  parameters  at  60  and  80  keVs,  compared  with 
the  exDerimental  data. 


A:  ■  data 


0  10  20  30  40  50  60  70  80  90 


Energy  (keV) 

Figure  18.  Simulation  of  transmission  spectra 
of  fiber  A,  fiber  B60  (B  with  length  of  60  mm) 
and  fiber  B30  (B  with  length  of  30  mm)  with 
their  best  fit  parameters  listed  in  Table  J, 
compared  with  the  experimental  data. 
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Figure  19.  Simulated  Transmission  spectra  of 
fiber  C,  fiber  D  and  fiber  E  with  their  best  fit 
parameters  listed  in  Table  I,  compared  with 
the  exnerimental  data. 
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1.  ABSTRACT 

Many  medical  imaging  and  industrial  applications  for  x  rays  require  large  area  optics  with  good 
scatter  rejection.  Preliminary  scatter  rejection  and  contrast  measurements  show  that  a  prototype  long 
borosilicate  optic  increases  the  contrast  by  a  factor  of  1 .7  by  decreasing  the  scatter  transmission  nearly  a 
factor  of  1 0  at  20  keV.  Since  borosilicate  optics  have  higher  scatter  transmissions  at  high  energies,  the 
optics  have  to  be  fairly  long  to  give  good  scatter  rejection  at  high  energies.  However,  long  optics  are 
complex  to  manufacture  and  have  increased  defect  rates.  Lead  glass  would  allow  the  optic  to  be  much 
shorter  and  still  give  good  contrast  enhancement,  because  of  the  superior  absorption  of  leaded  glass.  In 
order  to  investigate  the  feasibility  of  using  leaded  glass  polycapillary  optics  for  these  applications, 
measurements  and  simulations  have  been  performed  on  the  behavior  of  leaded  glass  polycapillary  fibers  in 
the  9-80  keV  energy  range.  These  fibers  have  an  outer  diameter  of  about  0.5  mm,  channel  diameter  of  12 
pm  and  fractional  open  area  from  40%  to  60%.  The  transmission  efficiencies  of  these  fibers  of  different 
types  and  lengths  were  measured  as  a  function  of  source  location  and  x-ray  energy.  The  measurements 
were  analyzed  using  a  geometrical  optics  simulation  program,  which  included  the  effects  of  roughness, 
waviness,  bending  and  glass  inclusions.  Despite  low  transmission  at  low  energies,  leaded  glass 
polycapillary  optics  with  a  length  of  30-60  mm  seem  promising  for  many  high  energy  (>20  keV)  x-ray 
applications.  Type  B  fibers  with  a  60  mm  length  have  a  transmission  efficiency  of  50%  in  the  35-40  keV 
energy  range,  which  is  83%  of  fiber’s  open  area,  and  very  low  scatter  transmission  of  0.1 8%>  at  80  keV.  The 
same  fiber  with  a  shorter  length  of  30  mm  has  a  transmission  which  is  higher  at  lower  energies,  and  more 
than  52%  in  the  20-40  keV  energy  range. 


2.  INTRODUCTION 

Polycapillary  optics,  consisting  of  arrays  of  tiny  hollow  glass  tubes  with  diameters  on  the  order  of  few 
microns,  is  a  relatively  new  technology  for  controlling  x-ray  beams.  A  polycapillary  fiber  is  a  capillary 
bundle  with  hundreds  or  thousands  of  hollow  channels,  as  shown  in  figure  1 .  X  rays  striking  the  interior  of 
the  glass  tubes  at  grazing  incidence  are  guided  down  the  tubes  by  total  external  reflection.  Arrays  of  curved 
or  tapered  capillaries  can  be  used  to  focus,  collimate  and  filter  x-ray  radiation.  Polycapillary  optics  have 

potential  usefulness  in  many  applications,  including  x-ray  lithography, 
medical  imaging,  crystallography  and  astronomy.1'3  Properties  of 
borosilicate  polycapillary  fibers,  fabricated  optics,  and  their 
applications  have  been  studied  extensively  in  the  8-80  keV  energy 
range.4'11  Borosilicate  single  fibers  with  a  length  of  140  mm  have  good 
transmission,  but  also  have  a  large  high-angle  scatter  transmission  of 
16%  at  energy  80  keV.12  Scatter  rejection  and  contrast  measurements 
using  a  tapered  borosilicate  glass  optic  are  presented.  Since  borosilicate 
optics  have  high  scatter  transmission  at  higher  energies,  the  optics 
would  need  to  be  very  long  to  provide  good  scatter  rejection  for  high- 
energy  applications.  However,  long  optics  lead  to  complexity  in 
manufacturing,  expense  and  increased  defect  rate.  For  many 
applications  of  scatter  rejection  optics,  large  areas  are  necessary  and 
large  area  long  optics  are  practically  difficult  to  produce.  Since  shorter 
optics  are  easier  to  manufacture,  optics  with  a  shorter  length  would 


Figure  1.  Photograph  of  a 
cross-section  of  a  leaded  glass 
polycapillary  fiber  with 
channel  diameter  of  12  pm  and 
fractional  open  area  of  60%. 
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provide  a  faster  route  to  bringing  the  benefit  of  polycapillary  x-ray  optics  to  these  applications.  Leaded 
glass  optics  are  advantageous  for  excellent  scatter  rejection  compared  to  borosilicate  glass  because  of  their 
superior  absorption.  To  explore  the  potential  of  lead  glass  scatter  rejection  optics,  the  transmission  for 
variety  of  leaded  glass  single  fibers  was  measured  in  the  energy  range  from  9  to  80  keV.  The  results  were 
analyzed  using  a  geometrical  optic  simulation  program.  The  use  of  the  simulation  to  analyze  the 
experimental  data  for  60  mm-!ong  type  B  fibers  is  demonstrated.  The  fitting  process  was  used  to  determine 
the  best-fit  parameters  of  the  simulation  for  five  other  leaded  glass  fibers. 


3.  BASIC  PRINCIPLE 


Capillary  x-ray  optics  use  multiple  total  external  reflections  to  guide  grazing  incidence  x  rays  in  arrays 
of  tiny  hollow  glass  tubes.  For  the  leaded  glass  capillaries  that  were  used  in  the  experiments  described  here, 
the  critical  angle  for  reflection  is13 

Gc  =  35  (mrad)  /  E  ( keV),  (1) 

The  higher  the  photon  energy,  the  smaller  is  the  critical  angle. 


4.  CONTRAST  ENHANCEMENT  USING  BOROSILICATE  POLYCAPILLARY 

OPTICS 


4.1  Experimental  Setup 


The  experimental  setup  for  scatter  rejection  and  contrast  measurements  is  shown  in  figure  2.  The  x-ray 
generator  used  in  the  experiment  was  a  Mo  source  with  a  300  pm  source  spot  size  and  maximum  operating 
voltage  of  35  kVp.  The  source  head  was  mounted  on  two  stages  which  could  be  translated  in  the  two 
directions  transverse  to  the  x-ray  beam.  Thick  lead  shielding  was  used  to  prevent  background  x  rays  from 
reaching  the  imaging  plate.  The  borosilicate  monolithic  optic  was  put  on  aluminum  plate  and  covered  by 
powder  filling  to  prevent  x-ray  leakage  around  the  optic.  A  Fuji  imaging  plate  with  50  pm  resolution  was 
used  to  record  the  image.  The  Lucite  phantom  used  in  the  contrast  measurement  was  45  mm  thick.  Several 
lead  strips  with  different  width  were  attached  to  the  phantom  to  measure  scatter  transmission. 

4.2  Optic 

A  prototype  borosilicate  monolithic  optic,  166  mm  in  length,  with  input  diameter  of  4.5  mm  and  7.6  mm 
output,  was  used  to  measure  scatter  fraction  and  contrast  enhancement  of  a  Lucite  phantom  with  a  thickness 
of  45  mm.  The  transmission  of  the  whole  optic  was  35%  at  20  keV  and  7%  at  60  keV.  A  central  part  of  the 
optic,  4  mm  in  diameter  at  the  output,  has  good  transmission  up  to  energy  40  keV. 
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Figure  3.  Output  uniformity  of  a  long 
monolithic  optic,  L=166  mm,  Dinpu,  = 
4.5  mm,  Doutput  =7.6  mm 


Figure  4.  Sketch  of  x-ray  paths 
for  scatter  measurement. 


The  output  uniformity 
of  the  optic  is  shown  in 
figure  3.  The  channels  of 
the  outer  part  of  the 
optic,  which  does  not 
transmit  well,  may  be 
blocked  by  glass 
inclusions.  Such  an  optic 
provides  magnification, 
and  thus  enhanced 
resolution,  as  well  as 
reducing  scatter. 


4.3  Scatter  fraction  Measurement 


To  measure  the  scatter  transmission  of  the  optic,  several  lead  strips  with  different  width  but  the  same 
thickness  were  attached  to  the  phantom,  as  shown  in  figure  4.  The  scatter  fraction,  SF,  was  calculated  using 
formula 


SF  =  Is  /  ( Is  +  Ip  ) 


(2) 


where  Is  is  the  intensity  behind  the  lead  strip  (scatter  intensity)  and  IP  is  the  primary  intensity,  which  is  the 
intensity  around  the  lead  strip  subtracted  by  scatter  intensity.  The  optic  reduced  the  scatter  fraction  by 
nearly  a  factor  of  10,  as  shown  in  figures  5a  and  5b. 


Figure  5a.  Scatter  fraction  with  optic.  Figure  5b.  Scatter  fraction  without  optic. 


4.4  Contrast  Measurements 


* 

i  depth,  d 

i 
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The  Lucite  phantom  used  in  the  contrast  measurements 
was  45  mm  thick  and  had  four  holes  with  the  same  diameter 
of  1.5  mm  and  different  depths,  from  10  mm  to  35  mm,  as 
listed  in  Table  I.  The  diameter  of  the  holes  were  smaller 
than  the  central  area  the  optic  which  transmits  well.  Each 
image  was  taken  on  a  single  hole.  A  sketch  of  the  Lucite 
phantom  with  two  of  the  holes  is  shown  in  figure  6.  The 
contrast  was  measured  with  and  without  optic.  The  images 
with  and  without  the  optic  are  shown  in  figures  7a  and  7b. 
The  contrast,  C,  is 


Figure  6.  Sketch  of  Lucite  phantom 


C  =  In  ( I'  / 1 ) 


(3) 
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where  I  is  the  intensity 
near  the  hole  and  I  is  the 
intensity  through  the  hole. 
The  contrast  with  and 
without  the  optic,  and  also 
contrast  enhancement,  are 
listed  in  Table  I.  The 
contrast  enhancement  is 

Figure  7b.  Lucite  phantom  the  ratio  of  the  contrast 
image  without  optic.  wlth  and  without  the  optic. 

The  results  show  that  the 

optic  increases  the  contrast  by  a  factor  of  1 .7  for  the 
shallowest  hole.  The  images  in  figures  7a  and  b  are  shown 
to  scale.  A  magnification  of  1 .8  was  obtained  with  the 
optic.  This  long  monolithic  optic  with  small  area  was  used 
to  magnify  the  image,  reduce  the  scatter  transmission  and 
enhance  the  contrast. 

4.5  Disadvantages  of  Long  Borosilicate  Optics 

Long  optics,  even  with  small  input  area,  tend  to 
have  defects  as  shown  in  figure  3.  A  long  optic  also 
leads  to  manufacturing  complexify  and  expense. 
Improvements  in  the  technology  of  fabricating 
monolithic  optics  have  been  achieved.14  Medical 
imaging  or  astronomy  applications,  however, 
require  large  area  optics.  For  high-energy 
applications,  borosilicate  optics  need  to  be  quite 
long.  Leaded  glass  would  allow  the  optic  to  be 
much  shorter  and  still  good  give  scatter  rejection 
because  of  strong  absorption  of  lead  glass.  The 
shorter  optics  are  easier  to  manufacture  and  could 
be  cheaper.  In  the  following  section,  several  leaded  glass  fibers  were  measured  and  analyzed  using  a 
geometric  simulation. 

5.  GEOMETRIC  SIMULATION 

In  the  geometrical  optic  simulation  used  in  this  article,  x-ray  transmission  through  hollow  glass  fibers  is 
simulated  by  tracing  a  large  number  of  x  rays  through  fibers.  The  simulation  includes  four  fiber  quality 
defects:  roughness;  waviness;  channel  blockage;  and  also  profile  corrections  that  can  be  approximated  by 
uniform  bending.  It  had  been  shown  previously  that  geometric  simulations  with  these  four  best  parameters 
are  in  good  agreement  with  experimental  data  for  borosilicate  single  fibers  at  1-  80keV.10'  15 

The  geometric  algorithm  is  a  two  dimensional  approximation.  Because  of  the  small  critical  angle,  the 
velocity  of  the  x-ray  photon  along  the  capillary  axis  is  nearly  equal  to  c,  the  velocity  of  light.  The  trajectory  of 
an  x  ray  can  be  reduced  to  the  two  dimensional  motion  inside  the  capillary  cross  section. 1(1  Deviations  of  the 
channel  from  a  straight  path  are  approximated  by  uniform  bending  described  by  a  bending  radius,  R.  The 
bending  of  the  channel  makes  the  apparent  motion  of  the  x  ray  within  the  channel  cross  section  similar  to  that 
of  a  classically  accelerated  particle. 

Surface  roughness  is  formally  parameterized  by  a  correlation  length  and  root-mean  square  displacement 
of  the  rough  surface,  z.  Because  changes  in  the  correlation  length  change  the  effective  roughness  height,  and 
due  to  the  lack  of  actual  correlation  data,  the  correlation  length  was  fixed  at  6  pm  for  all  simulations. 

Changing  this  value  would  cause  a  proportional  scaling  all  of  reported  roughness  heights.  The  value  was 
chosen  for  borosilicate  fibers  so  that  the  roughness  heights  agree  with  atomic  force  microscopy  data.10  The 
surface  roughness  decreases  the  apparent  reflectivity  of  the  channel  walls  and  therefore  the  transmission  in 
capillary  channels.  The  third  capillary  surface  quality  parameter  is  waviness,  which  occurs  on  a  spatial  range 
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FigurPba.  Lucite  phantom  image  with 
optic. 


Table  I.  Contrast  and  contrast  enhancement 
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shorter  than  the  capillary  length  and  longer  than  roughness.  The  average  effect  of  the  waviness  can  be 
considered  as  a  random  tilt  of  glass  wall.  In  the  simulation,  these  tilt  angles  are  normally  distributed  with  a 
mean  value  of  zero  and  a  standard  deviation  w.  The  fourth  polycapillary  defect  is  channel  wall  blockage, 
which  is  modeled  by  a  filter  layer  of  glass  of  thickness  t. 


6.  MEASUREMENTS  OF  LEADED  GLASS  FIBERS 

6.1  Experimental  Apparatus 

The  experimental  apparatus  is  shown  in  figure  8. 

6.1.1.  Source 

The  x-ray  generator  used  in  the  experiment  was  a  low  current  Microfocus  MS50  with  a  50  pm  spot  size, 
tungsten  target,  and  a  maximum  operating  voltage  of  100  kVp.  The  source  head  was  mounted  on  a  stage 
which  could  be  moved  in  the  two  directions  transverse  to  the  x-ray  beam. 

In  order  to  reduce  the  background  of  x  rays,  this  movable  source  was  enclosed  in  a  6-mm-thick  lead 

box  with  a  100x175  mm  aperture  in  the  front. 

The  source  to  fiber  distance  was  1 1 00  mm.  In  the 
measurement  for  photons  with  energy  higher 
than  50  keV,  a  5-mm-thick  aluminum  plate  was 
used  as  a  filter  between  the  source  and  the  fiber 
to  remove  the  low  energy  photons,  and  reduce 
the  dead  time  in  the  detector. 

6.1.2.  Pinhole 

A  200  pm  diameter  pinhole  through  4  mm  of 
lead,  2  mm  tungsten,  and  1  mm  of  tantalum 
(attached  together),  was  placed  5mm  away  from 
the  fiber  at  the  entrance  end.  The  pinhole  was 
attached  to  a  6-mm-thick  lead  shield.  The 
pinhole  and  lead  shield  together  block  scattered  x 
rays  and  leakage  around  the  outside  of  the  fiber. 
The  pinhole  was  smaller  than  the  fiber  but  still  covers  hundreds  or  more  channels.  Its  small  size  was  chosen 
not  only  to  avoid  the  leakage  around  the  fiber  and  reduce  dead  time  in  the  detector,  but  also  to  keep  the 
source-fiber  distance  as  small  as  possible  while  keeping  the  photon  entrance  angle  at  the  edges  of  the  pinhole 
to  less  than  the  critical  angle  at  energies  up  to  80  keV. 

6.1.3.  Fiber 

The  fiber  was  held  straight  by  a  finely  machined  groove  in  an  aluminum  plate  and  covered  by  iron 
powder  to  prevent  x-ray  leakage  around  the  fiber.  The  aluminum  plate  was  also  mounted  on  a  stage  which 
can  be  translated  in  two  orthogonal  directions  transverse  to  the  beam.  All  stages  were  mounted  on  rail 
carriers,  which  could  be  moved  along  the  beam  directions.  Five  kinds  of  fibers,  listed  in  Table  II,  were 
measured.  The  outer  cross  sections  of  all  the  fibers  are  roughly  hexagonal. 

6.1.4.  Detector 

The  detector  was  a  high  purity  germanium  detector  with  about  200  eV  resolution  at  5.9  keV  and  550  eV 
at  122  keV.  The  detector  was  placed  behind  the  fiber.  The  distance  between  the  source  and  detector  was 
fixed  to  keep  air  absorption  constant.  Motion  control  and  data  collection  were  all  controlled  with  a  small 
computer. 


X-ray  source 


Fiber  held  by 
metal  plate 


Lead  shielding  with  200 
pm  pinhole 


Detector 


Figure  8.  Experimental  setup  for  fiber 
measurement. 


7.  RESULTS  AND  ANALYSIS  FOR  LEADED  GLASS  FIBERS 


7.1.  High-angle  Transmission 
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Figure  9.  Normalized  intensity  vs  source  position, 
taken  by  moving  the  source  in  the  horizontal 
direction  perpendicular  to  the  fiber  with  fiber 
and  pinhole  fixed. 


O 


Figure  10.  Large  angle  (cut  through)  transmission 
versus  energy,  measured  with  the  source  moved 
away  from  the  aligned  position  to  an  angle  larger 
than  the  critical  angle.  Dash  lines  are  theoretical 
calculation  using  equation  4. 


The  pinhole  was  located  in  the  center  of  the  fiber  to  avoid  leakage.  Figure  9  shows  a  source  scan  plot  of 
transmission  as  a  function  of  source  location  for  eight  different  energy  windows  of  approximately  2  keV  in 
width.  The  transmission  plots  are  narrower  for  the  higher  energy  windows,  as  expected  due  to  the  decrease  in 
critical  angle  with  photon  energy.  For  leaded  glass  fibers,  because  of  the  high  absorption  of  lead,  the 
background  signal  from  photons  with  incidence  angle  greater  than  the  critical  angle  is  very  small.  Most  of 
these  photons  are  absorbed  rather  than  cutting  through  the  glass  wall.  The  cut  through  transmission  is  3.6%  at 
80  keV  for  a  30-mm-long  type  B  fiber  and  0.18  %  for  a  60-mm-long  fiber.  Measured  and  theoretical  high- 
angle  transmission  at  several  energies  are  shown  in  figure  10.  The  theoretical  values  are  calculated  as 


Tc  =  exp  {-L  (  1-f)  pp  }  (4) 

where  L  is  the  fiber  length,  f  is  the  fractional  open  area  of  the  fiber,  p  is  the  density  of  the  glass,  and  p  is  the 
mass  absorption  coefficient  of  the  glass  obtained  from  the  tabulated  values.16  The  transmission  results  which 
follow  are  channel  transmissions  from  which  the  very  small  high  angle  cut  through  and  fluorescence 
background  have  been  subtracted. 


7.2.  Simulation  Analysis 

Table  I.  Fiber  description  and  simulation  parameters 
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Figure  11.  Simulated  transmission  spectra  of  a 
type  B60  fiber  for  different  bending  radii  are 
compared  with  the  experimental  data.  The  figure 
shows  that  bending  alone  can  not  fit  the  data  since 
the  smallest  bend  underestimates  the  transmission 
at  high  energies. 


Figure  12.  Simulation  of  source  scan  curves  with 
different  bending  compared  with  the  experimental 
data  at  9  keV.  The  simulated  curve  does  not  become 
narrow  enough  even  for  a  bend  radius  of  R=20  m. 


Transmission  for  the  leaded  glass  fibers  in  Table  II  was  measured  in  the  energy  range  9-  80keV.  To 
illustrate  the  effect  of  each  simulation  parameter,  the  application  of  the  simulation  to  the  experimental  data 
for  fiber  type  B60  (a  60  mm-long  type  B  fiber)  is  discussed  in  this  section. 


7.2.1.  Bending  effects 

A  slight  bending  can  dramatically  reduce  the  transmission  at  the  high  energies  because  of  the  small 
critical  angles.  The  measured  transmission  as  a  function  of  energy  along  with  the  simulation  of  different 
radii  for  a  type  B60  fiber  is  shown  in  figure  11.  The  simulation  with  bending  alone  can  not  fit  the 
experimental  data  since  at  smallest  radius  the  transmission  is  too  high  at  lower  energies  and  too  low  at  high 
energies.  Since  any  waviness  correction  will  reduce  the  transmission,  the  radius  of  bend  has  to  be  larger 
than  32  m.  Simulations  with  different  bending  radii  and  the  experimental  transmission  source  scan  at  9  keV 
are  shown  in  figure  12.  Even  with  a  sharp  bend,  the  simulated  transmission  source  scan  at  9  keV  is  not 
narrow  enough. 

7.2.2.  Waviness  effect 

Simulations  with  different  waviness  corrections  are  compared  to  the  experimental  data  of  transmission 
spectra  in  figure  13.  The  waviness  alone  can  not  fit  the  data.  The  simulation  shows  that  the  waviness  causes 
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Figure  13.  Simulated  transmission  spectra  are 
compared  with  the  experimental  data.  The  waviness 
correction  causes  the  transmission  to  drop  primarily 
in  the  middle  energy  range. 


9  keV :  ■  exp 


Source  Position(mm) 
Figure  14.  Simulated  source  scan  curves  with 
different  waviness  are  compared  with  the 
experimental  data  at  9  keV.  A  waviness  of  0.5 
mrad,  which  is  too  large  at  higher  energies,  is 
needed  to  fit  the  data. 
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Figure  15.  Simulated  source  scan  curves  with  different 
waviness  are  compared  with  the  experimental  data  at  40 
keV.  Since  the  simulated  curve  with  a  waviness  of  0.15 
mrad  fits  the  data  fairly  well,  it  is  used  as  the  value  for 
waviness  in  the  best  fit  simulation  parameters  in  Table  I. 
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Figure  16.  Simulated  transmission  spectra  for  a 
type  B60  fiber  for  different  roughness  corrections 
are  compared  with  experimental  data. 


the  transmission  to  drop  rapidly  at  medium  energies  but  not  rapidly  enough  at  higher  energies.  Figure  13 
shows  that  the  waviness  correction  has  to  be  smaller  than  0.3  mrad  to  fit  the  data.  Simulated  transmission 
source  scans  with  different  waviness  are  compared  to  the  experimental  data  at  9  keV  in  figure  14.  The 
simulation  with  a  waviness  of  0.5  mrad  fits  the  spectrum  source  scan  fairly  well  at  9  keV.  However,  that 
was  seen  to  be  too  large  for  the  transmission  spectrum  of  figure  13.  Figure  15  compares  the  experimental 
data  to  simulations  at  40  keV,  and  shows  that  the  maximum  value  of  waviness  is  0.15  mrad. 


7.2.3  Roughness  effects 

Simulations  of  the  transmission  of  fiber  type  B60  with  and  without  roughness  corrections  are 
compared  to  the  experimental  data  in  figure  16.  For  a  perfect  capillary,  the  transmission  spectrum  is  flat, 
since  the  distance  between  the  source  and  the  fiber  is  large  enough  to  ensure  that  every  photon  which  hits  the 
glass  wall  has  grazing  angle  smaller  than  its  critical  angle.  The  transmission  is  almost  flat  at  all  energies,  but 
for  high  roughness  drops  more  rapidly  at  medium  energies  than  at  high  energies.  The  transmission  drop 
depends  on  both  the  surface  reflectivity,  which  is  reduced  by  roughness,  and  the  number  of  reflections  that 
the  photon  undergoes  through  the  channel.  The  transmission  is  proportional  to  R"  where  R  is  reflectivity,  and 
n  is  the  number  of  bounces.  For  a  straight  capillary,  the  average  number  of  bounces  of  the  photon  along  the 
channel  is  small,  usually  less  than  three.  Thus  the  roughness  has  only  a  fairly  small  effect  on  the 
transmission  spectrum  in  figure  16.  The  transmission  as  a  function  of  source  location  at  9  keV  is  shown  in 
figure  17.  Because  photons  experience  more  than  21  reflections  for  large  source  displacements,  the 
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Figure  17.  Simulated  source  scan  curves  for  a 
type  B60  fiber  for  different  roughness  are 
compared  with  experimental  data  at  9  keV. 
The  curve  is  too  low  at  large  displacements  for 
the  highest  value  of  roughness. 


Energy'  (keV) 


Figure  18.  Simulated  transmission  with  different 
layer  thicknesses  alone  compared  with  experimental 
data.  The  figure  shows  that  the  maximum  layer 
thickness  for  a  type  B60  fiber  is  35  pm  to  fit  the  low 
energy  transmission. 
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Figure  19.  Simulated  transmissions  with  and 
without  leaded  glass  filter  for  type  B60  fiber  are 
compared  with  the  experimental  data.  Best  fit  is 
simulated  transmission  with  R=56  m,  w=0.15  mrad, 
z=1.7  nm,  and  t=33  pm. 


Figure  20.  Simulation  of  source  scan  curves  with 
z=1.7  nm  ,  R=56  m,  and  w=0.15  mrad,  and  with 
roughness  z=1.7  nm  alone,  compared  with  the 
experimental  data.  Since  the  actual  transmission 
at  this  energy  is  quite  low,  the  source  scan  curve 
might  be  somewhat  distorted,  and  thus  the 
simulation  does  not  fit  the  data  well. 


transmission  at  high  angles  is  sensitive  to  roughness.  The  transmission  is  too  small  if  the  roughness  is  larger 
than  1.7  nm. 


7.2.4.  Channel  blockage 

By  modeling  channel  blockage  with  a  glass  fdter  layer  of  thickness  t=  35  pm,  the  simulated  transmission 
spectrum  with  the  glass  layer  alone  fits  the  data  at  low  energies.  The  simulated  transmission  spectra  with 
different  layer  thicknesses  are  compared  with  the  experimental  data  in  figure  18. 

To  fit  the  high  energy  transmission  spectrum  data,  finally,  we  combine  the  waviness  determined  from 
source  scan  at  40  keV,  the  roughness  from  source  scan  at  9  keV,  and  add  the  bending  radius,  that  must  be 
larger  than  32  m,  to  achieve  the  simulation  shown  in  figure  19.  To  fit  the  experimental  data  for  whole 
spectrum,  the  simulation  must  also  include  absorption  by  a  leaded  glass  layer,  which  models  glass  debris  in 
the  channels.  As  shown  in  figure  19,  the  inclusion  of  a  33  pm  thick  leaded  glass  layer  gives  a  good  fit  with 
the  experimental  transmission  spectrum. 

The  best  fit  simulation  compared  to  the  experimental  source  scan  curves  are  shown  in  figures  20-23. 
The  simulated  transmission  source  scan  curves  fit  the  data  fairly  well  at  energies  from  20  to  80  keV.  The 
source  scan  at  9  keV  does  not  fit  the  simulated  transmission,  as  shown  in  figure  20.  This  may  be  partly 
because  the  actual  value  of  the  transmission  at  9  keV,  as  shown  in  figure  19,  is  quite  low. 

7.2.5  Simulation  summary 

The  primary  effect  of  bending  is  to  reduce  the  transmission  at  high  energies.  Channel  blockage  reduces 
the  transmission  at  low  energies.  Roughness  reduces  the  width  of  source  scan  at  low  energies.  Waviness 
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Figure  21.  Simulated  scan  curves  compared  with 
the  experimental  data  at  three  photon  energies. 


Figure  22.The  simulation  of  scan  curves  with 
best  fit  parameters  at  50  and  70  keV,  compared 
with  the  experimental  data. 
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primarily  reduces  the  transmission  and  the  source 
scan  width  at  mid  range  energies. 


7.3.  Fiber  Results 

Using  the  fitting  processes  described  above,  the 
best  fit  simulation  parameters  for  fibers  A,  C,  B60, 
D  and  E,  which  have  length  60  mm,  and  also  fiber 
B30,  which  has  length  30  mm,  were  obtained.  The 
best  fit  simulation  parameters  are  listed  in  Table  II. 
The  best  fit  simulated  transmissions  along  with  the 
experimental  data  are  shown  in  figures  24  and  25. 
The  outer  diameter  of  the  fibers  are  in  the  0.51-  0.53 
mm  range.  These  fibers  are  thin  and  flexible,  and 
therefore  are  difficult  to  keep  straight  in  the 
measurement  apparatus,  resulting  in  sharp  bending 
radii  and  relatively  poor  high-energy  transmission. 

Fiber  type  B  has  high  open  area,  low  channel  blockage,  and  the  least  bending  and  waviness.  Thus  this 
fiber  has  better  transmission  than  other  fibers.  Figure  24  shows  that  the  transmission  of  this  fiber  type  with 
a  shorter  length  (B30)  is  higher  at  lower  energies  because  of  its  thinner  blockage  layer.  This  is  consistent 
with  a  statistically  random  model  of  glass  inclusion.  Because  of  its  high  flexibility,  the  shorter  fiber 
experiences  more  bending  than  the  longer  fiber,  and  thus  the  transmission  falls  off  faster  at  higher  energies. 
The  shorter  fiber  has  a  higher  scatter  cut  through,  as  shown  in  figure  10.  The  different  open  area,  bending, 
waviness  and  blocked  channels  cause  the  differences  in  the  transmission  among  the  fibers.  Fiber  type  C  and 
D  have  the  same  open  area,  but  fiber  type  C  has  thinner  glass  layer,  as  listed  in  Table  II.  As  a  result,  fiber 
type  C  has  higher  transmission  at  lower  energies  than  fiber  type  D,  as  shown  in  figure  25. 


c 


Figure  23.  The  simulation  of  scan  curves  with  best 
fit  parameters  at  60  and  80  keV,  compared  with  the 
experimental  data. 


8.  CONCLUSIONS 

The  borosilicate  optic  increased  the  contrast  in  a  Lucite  phantom  image  by  a  factor  1 .7  by  removing  the 
scatter.  However,  to  provide  good  scatter  rejection  at  high  energies,  borosilicate  optics  have  to  be  fairly 
long  which  leads  to  an  increased  defect  rate.  The  strong  absorption  of  leaded  glass  allows  the  optics  to  be 
much  shorter.  Measurements  were  performed  on  leaded  glass  fibers  to  asses  the  potential  for  leaded  glass 
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Figure  24.  Simulation  of  transmission  spectra  of 
fiber  A,  fiber  B60  (B  with  length  of  60  mm)  and 
fiber  B30  (B  with  length  of  30  mm)  with  their  best 
fit  parameters  listed  in  Table  II,  compared  with 
the  experimental  data. 
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Figure  25.  Simulated  Transmission  spectra  of 
fiber  C,  fiber  D  and  fiber  E  with  their  best 
parameters  listed  in  Table  II,  compared  with  the 
experimental  data. 
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scatter  rejection  optics.  The  high  measured  transmission  (>45%)  for  fiber  type  B  shows  that  the 
performance  of  this  leaded  glass  fiber  is  quite  good  in  the  25-50  keV  energy  range.  The  geometric 
simulation,  with  roughness,  waviness,  bending  and  channel  blockage  corrections,  performs  well  in 
explaining  the  measured  transmission  spectra  above  15  keV.  According  to  this  simulation,  bending  is  the 
most  harmful  at  higher  energies,  and  glass  inclusions  that  block  the  channels  are  harmful  at  lower  energies. 
For  type  B  fibers,  even  though  the  roughness  is  high,  the  waviness  is  low  and  the  bending  radius  is  the 
largest  among  the  fibers.  The  high-energy  transmission  performances  could  be  improved  by  making  the 
fiber  diameter  larger  and  therefore  more  rigid  to  minimize  unintentional  bending.  Low  energy  performance 
is  improved  by  using  shorter  fibers.  Because  of  the  good  absorption  of  lead,  the  high-energy  photons, 
which  are  not  reflected,  are  absorbed  completely  by  60  mm-long  fibers,  even  at  80  keV.  The  high-angle 
transmission  at  80  keV  is  only  0.2%.  This  high-angle  transmission  is  very  much  less  than  the  16%  for  a 
longer  (140  mm)  borosilicate  fiber.  Thus  short  leaded  glass  polycapillary  x-ray  optics  are  very  promising 
for  scatter  rejection  in  high-energy  applications. 
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